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ABSTRACT 


This  study  was  conducted  in  a  natural  area  which  was  located 
approximately  40  miles  east  of  Edmonton,  Alberta.  The  natural  area  was 
an  island  that  was  situated  at  the  northern  end  of  Beaverhill  Lake. 

The  study  was  initiated  to  determine  the  herbage  production 
of  the  area  following  the  invasion  of  woody  species.  To  test  the 
significant  changes  of  herbage  production,  the  community  structure  of 
the  island  vegetation  was  determined. 

The  vegetation  on  the  island  was  sampled  during  1970  using 
conscious  stand  selection.  However,  it  was  found  that  the  successional 
nature  of  the  island  vegetation  rendered  conscious  stand  selection 
unsuitable.  In  1971,  a  restricted  random  method  of  plot  location  was 
employed  to  obtain  a  more  representative  sample  of  the  vegetation. 

Three  study  units  were  delineated  for  sampling.  These  were  the  recon¬ 
naissance,  gradient  and  island  study  units. 

In  this  study,  the  ordination  and  cluster  analysis  techniques 
have  been  combined.  The  vegetation  analysis  technique  resulting  from 
this  combination  has  been  termed  orainational  classification. 

The  analysis  of  these  study  units  revealed  that  clusters  of 
vegetation  were  easily  defined  and  the  relationships  among  these 
clusters  were  clearly  presented  in  orthogonal  cluster  ordination.  It 
was  found  that  the  cluster  distribution  along  the  first  (Z)  orthogonal 
axis  was  significantly  correlated  with  elevation,  soil  moisture,  tree 
and  shrub  density,  depth  of  organic  soil  and  the  amount  of  sodium  in  the 
soil.  The  second  (Y)  orthogonal  axis  was  significantly  correlated  with 
June  soil  temperature,  forb  weight  and  annual  herbage  weight  while  the 
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third  (Z)  orthogonal  axis  was  significantly  correlated  with  amount  of 
available  phosphorous,  elevation,  soil  pH  and  soil  moisture. 


The  X-axis  defined  the  clusters  which  fell  into  the  upland 
and  lowland  zones  while  the  Y  and  Z  axes  separated  the  clusters  within 
the  upland  and  lowland  zones,  respectively.  Three  vegetation  associations 
were  formed  from  the  ordinations  and  field  data.  These  were: 

(a)  the  Cirs-ium  -  Soivpus  shoreline  association, 

(b)  the  Calamagvost'is  -  Carex  low-moor  association,  and 

(c)  the  PopuZus  -•  SaZix  woodland  association. 

A  fourth  association,  the  Poa  -  MeZZZctiis  grassland  association,  was 
formed  when  the  vegetation  was  intensively  sampled. 

It  was  found  that  the  annual  herbage  production  of  the  PopuZus  ~ 
SaZix  woodland  association  was  significantly  less  than  the  other  assoc¬ 
iations,  The  invasion  of  woody  species  reduced  the  annual  herbage 
production  by  68%.  The  highest  forb  weight  was  produced  in  the  Cirsium  ~ 
Soivpus  shoreline  association  while  the  Co, Zamagvos vis  -  Capex  association 
had  the  highest  grass  and  grass-like  weight.  Litter  accounted  for  62% 
of  the  available  forage  on  the  island. 
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INTRODUCTION 


The  scientific  study  of  rangelands  in  Alberta  has  long  been 
associated  with  the  short  and  mid-grass  prairie  regions  (Wroe,  1971). 
Recent  works  by  Hilton  (1970)  and  Wroe  (1971)  have  been  conducted  in 
the  Alberta  parkland.  Moss  (1955)  described  the  parkland  as  Ma  mosaic 
of  prairie  patches  and  aspen  groves  with  prairie  occupying  the  drier 
situations  and  aspen  the  more  moist  and  sheltered  places."  It  is  the 
transition  zone  between  the  boreal  forests  to  the  north  and  west  and 
the  mixed  prairie  to  the  southeast.  Pettapie.ee  (1969)  reported  that 
this  belt  may  range  from  stands  of  almost  pure  Populus  tremuto'td&s^' 
forests  on  the  northern  edge  to  large  areas  of  prairie  with  a  few 
trees  around  the  sloughs  near  the  southern  boundary.  Bird  (1961) 
mentioned  that  the  northern  portion  was  dominated  by  Populus  tvenru ~ 
loides  while  the  southern  half  by  Festv.oa  scdbvella  grasslands. 

The  search  for  agricultural  land  has  left  few  areas  of  the 
parkland  of  Alberta  unharvested.  Natural  undisturbed  vegetation  has 
become  a  rarity  and  the  remaining  stretches  of  virgin  land  should  be 
recorded  and  preserved  as  a  living  record  of  our  botanical  past.  The 
study  of  these  natural  areas  attempts  to  provide  guidelines  for  the 
proper  management  of  parkland  vegetation,  preserving  its  character 
for  future  generations  while  simultaneously  harvesting  its  resources 
for  food  production. 

The  objectives  of  this  study  were: 

(a)  to  investigate  the  herbage  production  of  an  undisturbed 
portion  of  the  parkland  belt, 

(b)  to  obtain  natural  ecological  groups  of  vegetation,  and 

_ 

J'The  nomenclature  of  the  plants  follows  Moss  (1959)  and  the  common 
names  have  been  given  in  Appendix  5« 
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(c)  to  statistically  analyse  differences  among  these  veget- 
ational  groups  for  herbage  production  and  several  environmental  and 
edaphic  factors. 
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DESCRIPTION  OF  THE  AREA 


I.  Location  and  Size 

The  general  study  area  which  surrounded  Beaverhill  Lake  was 
approximately  40  miles  east  of  Edmonton,  Alberta.  The  receding  water 
level  of  Beaverhill  Lake  has  exposed  approximately  10,000  hectares  of 
land.  The  Government  of  Alberta  has  maintained  this  land  for  the  crown 
and  has  been  leasing  it  to  the  local  farmers  for  grazing  purposes.  The 
exposed  lake  bed  has  been  subjected  to  intensive  grazing  with  the 
exception  of  two  islands  on  the  northern  end  of  the  lake.  The  smaller 
island,  10  hectares  in  size,  has  been  the  nesting  colony  of  the  Ringbill 
gull  {Lavas  delawarensis  Ord)  and  the  California  gull  ( Lotus  ca lifomi cus 
Lawrence) .  The  larger  island  has  not  been  recently  disturbed  and  was 
selected  as  the  specific  study  area.  This  island,  approximately  50 
hectares  in  size,  was  approximately  1110  meters  long  and  450  meters 
wide.  The  legal  description  was  the  north  half  of  section  24, 
township  52,  range  18,  west  of  the  fourth  meridian.  An  aerial  photograph 
of  the  island  has  been  presented  in  Figure  1. 

II.  History 

The  Water  Resources  Division  of  the  Alberta  Department  of 
Agriculture  reported  the  present  level  of  Beaverhill  Lake  to  be  664.2 
meters  above  sea  level.  The  maximum  level  was  reported  in  1915  as  667.3 
meters  (Anonymous,  1915).  During  the  past  80  years  there  appeared  to 
have  been  a  cyclic  pattern  of  low  and  high  water  levels.  Early  residents 
of  the  area  reported  that  the  lake  was  nearly  dry  during  the  early  1890 ’s 
(M.  Eleniak,  personal  communication).  By  1915,  the  lake  had  returned  to 
its  maximum  level  and  then  again  receded  to  a  low  level  of  approximately 
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Figure  1*  Aerial  photograph  of  the  study  area. 
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663.0  meters  by  1949.  During  this  latter  dry  period,  gravel  was  hauled 
from  the  eastern  end  of  the  study  island  and  hay  was  harvested  around 
the  island  (M.  Eleniak,  -personal  communication) . 

III.  Geology  and  Soils 

Bowser  et  al.  (1962)  reported  that  Beaverhill  Lake  was  located 
above  the  Bearpaw  formation  with  the  northeastern  edge  of  the  lake  over 
the  Pale  Beds  formation.  The  predominant  parent  material  around  the  lake 
was  glacial  till  while  small  pockets  of  fine  textured  lacustrine  material 
and  medium  textured  alluvial  lacustrine  material  were  found  along  the 
northwestern  end.  The  predominant  soil  surrounding  Beaverhill  Lake 
was  the  Black  Solodized  Solonetz  subgroup.  The  exposed  lake  bed  was 
composed  of  coarse  textured  outvash  material  and  was  mapped  as  being  in 
the  Regosolic  Order.  The  beach  sand  soils  of  the  lake  bed  were  considered 
non-arable  and  were  not  described  by  Bowser  et  al.  (1962). 

IV.  Climate 

Bowser  et  al.  (1962)  reported  that  the  climate  of  the  area  was 
continental  with  relatively  warm  summers  and  cold  winters.  The  May  to 
September  mean  temperature  was  56°F  while  the  mean  winter  temperature, 
November  to  March,  was  16°F.  July  was  the  warmest  month  with  a  mean 
temperature  of  61 °F  while  January  was  the  coldest  with  a  mean  temperature 
of  6°F.  The  frost-free  period  was  approximately  100  days  in  length. 

The  mean  annual  precipitation  was  40  cm  of  which  70%  fell  as 
rain.  June,  July  and  August  were  the  months  of  highest  rainfall  accounting 
for  one-naif  of  the  total.  The  average  snowfall  in  the  area  was  about 
127  cm.  The  November  to  March  precipitation  was  almost  exclusively  snow. 
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The  predominant  wind  direction  was  from  the  northwest  with  an 
average  velocity  of  16  kilometers  per  hour. 


LITERATURE  REVIEW 


I.  Marsh  Vegetation 

Moss  (1955)  reported  that  the  mesic  to  xeric  sequence  on  the 
borders  of  alkaline  sloughs  in  the  Alberta  parkland  was  Typha  -  Scirpus 
to  Scolochloa  -  Carex  atherod.es  to  Glyceria  -  Carex  spp.  to  Populus 
balsami-fera .  Lewis  et  al.  (1928)  studying  reed-swamps  in  Central 
Alberta,  found  that  the  common  succession  for  non-saline  large  sloughs 
was  Typha  -  Scirpus  to  Scolochloa  -  Carex  atherodes  to  Glyceria  -  Carex 
spp.  to  Calamagrostis  americana  -  Poa  palustris  to  Calacmagrostis  cana¬ 
densis  -  Sa'iix  to  Populus .  They  distinguished  the  Typha  -  Scirpus  to 
Scolochloa.  -  Carex  atherodes  as  reed-swamp  while  the  remainder  excluding 
Populus  was  termed  low-moor.  Lewis  et  al.  (1928)  reported  that  the 
Scirpus  -  Typha  association  was  open  and  standing  in  several  inches  of 
water.  The  Scolochloa  -  Carex  atherodes  association  was  a  closed  commun¬ 
ity  surrounding  the  Scirpus  -  Typha  association.  They  also  noted  that 
when  the  low-moor  association  was  continuously  mowed  Juncus  balticus 
(ater)  and  Poa  pvatensis  were  common  and  that  the  Calamagrostis  amer- 
icaivi  association  occurred  on  the  drier  sites  of  the  low-moor  association. 
Moss  (1953)  found  Scolochloa  festucacea  and  Carex  atherodes  as  dominants 
of  reed-swamps  of  northwestern  Alberta  and  remnants  of  a  Calamagrostis 
inexpansa  -  Carex  atherodes  marsh. 

Walker  and  Coupland  (1968)  found  that  Scolochloa  festucacea 
showed  a  preference  for  slightly  saline  areas,  while  Carex  atherodes 
favored  non -saline  areas.  They  found  that  mechanical  disturbance  at  the 
edge  of  a  slough  resulted  in  the  increase  of  Eordeum  jubatum _,  Rumex  spp., 

i 

Sonchus  arvensis ,  Cirsium  sppf ,  Taraxacum  officinale s  Polygonum  ccoiculare 
and  Artemisia  biennis.  Grazing  favored  Poa  palustris ,  Glyceria  grandis 
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and  Bechmannia  syzigachne ;  when  moist,  mowing  favored  ScoZochZoa 
festucacea  and  Cave. x  athevodes ;  while  under  drier  conditions,  Poa 
paZustvis  was  favored. 

Bird  (.1961)  reported  that  common  vegetation  around  permanent 
or  semi-permanent  sloughs  was  Cavex  -  ScoZochZoa 3  Typha  and  Seivpus 
zones.  Fresh  water  sloughs  were  generally  surrounded  by  a  zone  of 
SaZix  spp,  with  PopuZus  tvemuZoides  surrounding  the  SaZix  spp. 

Dix  and  Smeins  (1967)  found  in  North  Dakota  that  Seivpus  spp. 
and  Typha  spp,  were  the  dominant  emergent  marsh  species  while  ScoZoehZoa 
festucacea  ana  Cavex  athevodes  were  commonly  dominant  where  standing 
water  was  present  for  long  periods.  The  authors  found  that  CaZamagvostis 
inexpansa  and  Cavex  Zasiocavpa  (Zanuginosa)  were  associated  in  depressions 
which  were  inundated  with  water  for  short  periods.  They  defined  marsh 
as  stands  of  Phvagmites  communis 3  Typha  ZatifoZias  and  Seivpus  spp. 
while  meadows  ware  dominated  by  ScoZocfoloa  festucacea 3  Cavex  athevodes 
and  PhaZcvis  avurAinaeea . 

II.  Mathematical  Analysis  of  Vegetational  Relationships 

The  structure  of  vegetation  has  been  largely  determined  by 
the  differential  response  of  species  to  a  set  of  environmental  factors 
(Beals,  1960j  Ayyad  and  Dix,  1964).  The  greater  the  similarity  in  the 
response  of  two  species  to  the  environmental  factors  of  any  site,  the 
greater  the  chance  the  two  species  will  have  of  being  present  at  that 
site.  Therefore  the  resemblance  among  sites  can  be  determined  by  their 
species  composition  (Swan  et  al.,  1969).  The  degree  of  resemblance 
among  sites  may  be  measured  using  various  coefficients  (Sokal  and  Sneath, 
1963) .  The  coefficients  used  in  plant  ecology  have  been  the  coefficients 
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of  association  and  of  distance. 

A.  Measures  of  Site  Resemblance 

The  most  commonly  used  coefficient  has  been  that  of  association 
(Cole,  1949;  1957;  Bray  and  Curtis,  1957;  Williams  and  Lambert,  1959; 

1961;  Beals,  1960).  Bray  and  Curtis  (1957)  and  Beals  (1960)  used  the 
Index  of  Similarity  (S)  to  derive  a  measure  of  the  association  between 
vegetational  sites.  The  Index  of  Similarity  formula  used  by  Bray  and 
Curtis  (1957)  and  Beals  (1960)  was  S  =  2W/ (A+B) ,  where  A  was  the  sum  of 
the  quantitative  species  values  (canopy  coverage,  frequency,  density, 
etc.)  for  one  site,  B  was  the  equivalent  of  the  other  site,  and  W  was 
the  sum  of  the  lesser  values  of  those  species  which  were  present  in 
both  sites.  These  authors  transformed  the  coefficient  of  association 
into  a  coefficient  of  distance  which  they  called  the  Index  of  Dissim¬ 
ilarity  (D)  by  subtracting  the  Index  of  Similarity  from  the  maximum 
Index  of  Similarity  (D  =  S  -  S)  found  in  the  study. 

Orloci  (1966)  presented  an  alternative  association  coefficient 
which  he  termed  the  weighted  similarity  coefficient.  This  was  the 
coefficient  of  principal  components  ”f actor”  analysis. 

The  coefficient  of  distance  has  been  most  extensively  used 
in  numerical  taxonomy  (Sokal  and  Sneath,  1963).  One  coefficient,  the 
Euclidean  distance,  has  recently  been  used  in  plant  ecology  (Austin 
and  Orloci,  1966;  Orloci,  1966;  Swan  et  al.,  1969).  The  coefficients 
of  distance  have  been  based  upon  a  geometric  model  (Euclidean  space) 
which  has  been  defined  by  a  set  of  orthogonal  axes  (Sokal  and  Sneath, 

1963).  In  plant  ecology,  individual  sites  vrere  located  according  to 
their  species  composition  in  a  Euclidean  space  of  dimension  (number  of 
orthogonal  axes)  equal  to  the  number  of  species  present  (Swan  et  al.,  1969). 
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The  spatial  position  of  individual  sites  relative  to  one  another  can 

be  obtained  by  an  estimate  of  the  intersite  distance  (Orloci,  1966). 

Orloci  (1966)  suggested  that  the  intersite  distance  was  obtained  using 

the  Pythagorean  theorem.  If  A  and  B  were  two  sites  with  axes  values 

(X^,  X0^, . . . ,  X,^)  and  (X^.,3,  X1>T. , . .  „ ,  X^g)  ,  respectively  for  N  orthogonal 

axes,  then  their  intersite  distance  (D^,)  was  computed  as: 

N 


-  r  Z  (X.  -  X,„)2]1/2. 


DAB  L '"iA  “IB 
1=1 


Since  each  orthogonal  axis  was  defined  by  quantitative  values  of  a 
species,  the  intersite  distance  may  be  computed  directly  from  the  value 
defining  the  site  position  on  each  axis  (Orloci,  1966).  Sokal  and 
Sneath  (1963)  reported  that  the  intersite  distance  increased  with  the 
number  of  species  used  and  suggested  the  following  correction: 


dab  ■  !<V2'Nll/2- 

This  has  been  defined  by  Sokal  and  Sneath  (1963)  as  the  average  distance. 
Another  distance  coefficient  used  in  numerical  taxonomy  was 
the  arithmetic  distance  (Carmichael  and  Sneath,  1969).  The  formula 
given  was: 


N 


IX  - 
1  xA 


X±B 


1=1 

which  had  the  same  notation  as  the  Euclidean  distance.  Sokal  and  Sneath 
(1963)  corrected  the  arithmetic  distance  by  dividing  it  by  the  number 
of  species  and  called  this  distance  the  mean  character  difference. 

They  suggested  that  this  distance  was  simpler  to  use  than  the  Euclidean 
distance  but  the  mean  character  difference  will  always  underestimate 
the  true  Euclidean  distance  between  sites  in  Euclidean  space  especially 
when  species  differences  were  small  while  others  were  large.  No  record 


in 


.  *  . 


•  .  ■> 


• 

• 

.• .  -  ; 

4* 

■  * 

*  ■  * 

•  '• 

* .  - 

'  '•  V-  -'-*  -  '  •  p 

•»  i3 

■U  - 

% 

r 

m 

- 

.1  * 

*  . 

•  V 

*  ;  .  ■ 

*  V* » 


-11- 


of  the  usage  of  the  arithmetic  distance  or  the  mean  character  difference 
could  be  found  in  plant  ecology  literature. 

When  all  possible  intersite  distances  were  computed  using 
any  one  of  the  aforementioned  coefficients,  a  symmetrical  matrix  was 
produced  in  which  the  main  diagonal  values  were  zeros  and  the  rows  and 
columns  were  defined  by  the  number  of  sites  compared.  This  intersite 
distance  matrix  was  the  basis  of  the  ordination  technique  and  some 
numerical  classification  techniques. 

B.  Ordination 

Two  continuous  structuring  techniques  used  in  plant  ecology 
have  been  principal  components  "factor”  analysis  (Goodall,  1954;  Austin 
and  Orloci,  1966)  and  ordination  (Bray  and  Curtis,  1957;  Beals,  1960; 
Gittens,  1965a;  1965b;  Austin  and  Orloci,  1966;  Orloci,  1966;  Swan  et  al., 
1969).  Ordination  as  defined  by  Swan  et  al.  (1969)  was  a  construction 
of  a  geometrical  model  which  was  analogous  to,  yet  simpler  to  construct 
and  comprehend  than  principal  components  "factor”  analysis.  Ordination 
has  been  the  most  popular  of  the  two  continuous  structuring  techniques 
in  plant  ecology. 

Lambert  and  Dale  (1964)  stated  that  ordination  was  a  technique 
whereby  each  site,  was  placed  in  relation  to  all  others  on  one  or  more 
axes.  Swan  et  al.  (1969)  assumed  that  the  distribution  of  many  species 
was  related  to  a  limited  number  of  environmental  factors  and  that  the 
axes  defined  during  ordination  should  reflect  these  factors.  This 
concept  was  also  suggested  by  Bray  and  Curtis  (1957),  Beals  (1960)  and 
Ayyad  and  Dix  (1964) ,  The  construction  of  the  ordination  axes  as  defined 
by  Beals  (1960),  Orloci  (1966)  and  Swan  et  al.  (1969)  were  all  basically 
similar;  however,  the  axes  construction  method  of  Swan  et  al.  (1969) 
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presented  the  "best  stand  [site]-def ined  axes.” 

Swan  et  al.  (1969)  chose  the  first  orthogonal  axis  as  the  one 
which  accounted  for  the  highest  proportion  of  the  sums  of  squares  of 
all  intersite  distances  derived  from  the  Euclidean  distance.  The  second 
axis  accounted  for  the  next  highest  proportion  of  the  intersite  sum  of 
squares  not  accounted  for  by  the  first  axis.  The  second  axis  was  per¬ 
pendicular  to  the  first.  Other  axes  were  constructed  from  the  residual 
sums  of  squares.  The  criterion  of  axis  selection  has  been  fully  dis¬ 
cussed  by  Swan  et  al.  (1969).  They  believed  that  their  method  was  simpler 
to  understand  than  the  analogous  but  more  sophisticated  principal 
components  "factor"  analysis  because  it  was  based  on  axes  derived  from 
simple  geometric  constructions.  Their  method  was  also  found  to  be 
suited  to  species  ordination.- 
C.  Numerical  Classification 

The  ordination  technique  arranged  the  sites  on  axes  with 
the  position  of  each  site  being  determined  by  its  species  composition. 

In  classification,  the  sites  were  arranged  into  groups  with  each  group 
having  a  common  species  composition  (Lambert  and  Dale,  1964).  Greg- 
Smith  (1964)  defined  classification  as  the  "arranging  of  stands  [sites] 
into  classes,  the  members  of  each  of  which  have  in  common  a  number  of 
characteristics  setting  them  apart  from  the  members  of  other  classes," 

In  plant  ecology,  two  methods  of  finding  vegetational  groups  from  a  matrix 
of  observations  have  been  presented;  a  dendogram  (Williams  and  Lambert, 
1959;  1961),  and  cluster  analysis  (Stringer,  1969).  Both  of  these 
methods  have  been  used  extensively  in  numerical  taxonomy  (Sokal  and 
Sneath,  1963;  Carmichael  et  al.,  1968;  Carmichael  and  Sneath,  1969). 

The  dendogram  classification  as  presented  by  Williams  and 
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Lam'bert  (1959,  1961)  was  a  qualitative,  monothetic,  divisive  and  hier¬ 
archical  classification  system  which  was  derived  from  either  a  normal 
or  an  inverse  analysis.  The  major  discontinuities  between  groups  of 
sites  were  identified  during  normal  analysis.  This  was  accomplished 
by  using  the  correlation  of  the  presence-absence  records  of  all  possible 
species  pairs.  Inverse  analysis  detected  discontinuities  in  groupings 
of  species  from  the  correlation  of  site  pairs.  It  was  found  by  'Williams 
and  Lambert  (1959)  that  the  site  groupings  could  actually  be  mapped  and 
were  analogous  to  habitat  discontinuities.  Species  groupings  were 
found  to  be  ,lmerely  abstractions  of  ecologically-related  species,  held 
together  by  their  defining  quadrat,  but  rarely  representing  the  total 
flora  of  a  particular  site  except  in  the  case  of  the  first  hierarchical 
division”  (Williams  and  Lambert,  1961).  Normal  analysis  was  a  measure 
of  differences  in  species  richness  while  inverse  analysis  was  a  measure 
of  differences  in  species  abundance  (Williams  and  Lambert,  1961). 

Carmichael  and  Sneath  (1969)  suggested  that  the  dendogram 
could  not  illustrate  the  relationships  between  non-adjacent  groupings. 
They  presented  a  quantitative,  polythetic,  agglomerative  and  semi- 
hierarchical  taxometric  classification  technique  which  presented  the 
relations  between  blusters  located  in  Euclidean  space.  Carmichael’s 
(1970)  classification  was  a  method ’where  agglomeration  was  renewed  at 
each  level  of  classification,  thus  partitions  formed  at  one  level  did 
not  restrict  those  formed  at  another  level." 

The  cluster  analysis  developed  by  Carmichael  et  al.  (1968) 
began  with  a  set  of  points  in  Euclidean  space  as  defined  by  a  set  of 
orthogonal  axes.  They  were  interested  in  determining  if  "natural 
clusters"  (groups  of  sites)  existed  within  a  Euclidean  space  by  deter- 
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mining  the  relateclness  of  the  sites.  They  concluded  that  sites  forming 
a  cluster  must  be  more  related  than  sites  in  differing  clusters  and  that 
the  relatedness  of  sites  could  be  determined  from  the  intersite  distance. 
Sites  falling  within  a  cluster  had  relatively  short  intersite  distances 
among  themselves,  and  relatively  long  intersite  distances  to  sites  in 
other  clusters.  Carmichael  and  Sneath  (1969)  displayed  their  classif¬ 
ication  in  a  form  of  circles  representing  the  clusters  and  lines 
representing  the  cluster  relations.  The  authors  were  able  to  portray 
the  cluster  relations  inherent  within  a  Euclidean  space  and  thus  defined 
a  classification  technique  which  complemented  ordination. 

Stringer  (1969)  used  ordination  and  cluster  analysis  in  a 
study  of  grasslands  in  Bauff,  Jasper  and  Water ton  Lakes  National  Parks 
in  Alberta.  He  used  the  Bray  and  Curtis  (1957)  model  for  ordination 
and  the  George  and  Carmichael  (1966)  cluster  analysis.  Ha  found  that 
the  cluster  analysis  verified  the  discontinuity  between  the  Waterton 
and  the  Banf f-Jasper  stands  as  indicated  by  three-dimensional  ordination. 
Stringer  (1969)  also  found  that  cluster  analysis,  although  complimentary 
to  ordination,  delimited  mors  accurately  the  centrally  located  stands. 

D.  Problems  in  Distance  Measurement 

Swan  et  al.  (1969)  encountered  problems  with  Euclidean  distance 
measurements.  They  found  that  "the  calculated  distance  between  two 
stands  increases  with  the  number  of  species  for  which  quantitative 
measures  are  available."  The  intersite  distance  between  the  end  sites 
may  be  less  than  the  distance  to  a  site  having  more  species  near  the 
center  of  the  gradient,  thus  distorting  the  vegetational  patterns  and 
obscuring  vegetation-environment  relationships.  This  distortion  may 
be  overcome  by  using  the  average  distance.  They  found  that  species 
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values  of  zero  did  not  increase  the  distance  between  sites  even  though 
species  populations  were  changing  and  that  sites  which  were  quantitatively 
and  floristically  alike  may  be  adjacent  or  separated  by  other  sites  along 
a  vegetational  gradient. 
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METHODS 


The  vegetation  which  was  delineated  for  sampling  has  been 
defined  in  the  present  study  as  a  study  unit .  Three  study  units, 
reconnaissance,  gradient  and  island,  were  selected  to  sample  the 
vegetation  of  the  island  using  various  methods  of  plot  location. 

I.  Reconnaissance  Study  Unit,  1970 
Seventy  line  transects,  five  meters  in  length,  were  located 
within  stands  of  five  representative  communities  which  were  selected 
following  reconnaissance  (Daubenmire,  1952)  of  that  portion  of  the 
island  illustrated  in  Figure  2.  The  canopy  coverage  (%)  class  midpoint 
following  Daubenmire  (1959)  was  recorded  for  all  species  occupying  a 
percentage  of  the  line  transect.  The  coverage  classes  used  were: 

Midpoint  Range 


3  = 

0  - 

5% 

15  = 

5  - 

25% 

37  - 

25  - 

50% 

63  = 

50  - 

75% 

85  = 

75  - 

95% 

98  = 

95  - 

100% 

Two  20x50  cm  plots,  one  randomly  located  in  each  2.5  m 
segment,  were  harvested  during  early  August  1970.  The  green  herbage 
was  harvested  from  one  plot  and  sorted  to  grass  and  grass-like  species 
and  to  forbs.  Total  herbage  was  harvested  from  the  other  as  a  measure 
of  the  standing  crop.  The  collected  herbage  was  air-dried  in  the  labor¬ 
atory  and  the  weights  of  grass  and  grass-likes,  forbs  and  standing  crop 
were  recorded  for  each  appropriate  plot.  The  annual  herbage  weight  was 
determined  from  the  sum  of  the  grass  and  grass-like  weight  and  the  forb 
weight.  One  year  later,  a  20x50  cm  plot  immediately  adjacent  to  the 
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Figure  2.  Location  of  the  reconnaissance  plots  on  the  island. 
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first  plot  was  harvested  in  early  August,  sorted  and  air-dried  as  above, 
and  the  grass  and  grass-like  weight,  forb  weight  and  annual  herbage 
weight  were  recorded.  All  weights  were  converted  to  kilograms  per  hec¬ 
tare  and  comparisons  between  1970  and  1971  herbage  production  were  made. 

II.  Island  and  Gradient  Study  Units,  1971 
A.  Plot  Location 

Daubenmire  (1968)  and  Greg-Smith  (1964)  suggested  that  con¬ 
scious  stand  selection  should  be  practised  to  efficiently  sample  veg¬ 
etation,  Since,  the  vegetation  on  the  island  had  been  rapidly  changing 
during  the  last  40  years,  stands  of  vegetation  were  difficult  to  delineate 
during  the  reconnaissance  because  relic  species  were  common.  It  was 
decided  In  1971  that  a  more  representative  sample  of  the  vegetation 
would  be  obtained  using  a  restricted  random  method  of  plot  location. 

1.  Island  study  unit. 

An  arbitrary  base-line  was  drawn  east  to  west  along 
the  north  shore  of  the  island.  Perpendicular  cut-line  transects  were 
randomly  located  within  50  rrf  intervals  along  the  base-line  beginning 
from  the  eastern  end.  Each  transect  was  divided  into  50  m  segments  and 
a  point  was  randomly  selected  within  each  segment.  The  plot  center  was 
located  seven  meters  east  or  west  of  this  point  if  the  transect  was 
located  greater  than  or  less  than  25  m,  respectively  within  the  base-line 
segment.  A  distance  of  seven  meters  was  chosen  to  avoid  the  cut-line 
effect.  Each  plot  was  a  circle  of  two  meters  in  diameter.  The  plots 
distributed  along  23  transects  formed  the  island  study  unit.  The 
distribution  of  these  plots  has  been  displayed  in  Figure  3. 


All  measurements  were  to  the  nearest  meter. 
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Figure  3.  .Plot  distribution  in  the  island  study  unit 
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2.  Gradient  study  unit. 

Since  the  task  of  collecting  environmental  and 
edaphic  data  from  185  plots  would  have  been  astronomical,  two  transects 
(numbers  8  and  10)  were  selected  for  an  intensive  sample  of  the  major 
vegetational  gradient.  The  plots  were  randomly  selected  within  10  m 
segments  along  each  transect  and  were  located  seven  meters  west  and  east 
of  transect  numbers  8  and  10,  respectively.  The  island  study  plots  were 
located  first, while  the  gradient  study  plots  were  superimposed.  Data 
collected  for  an  island  study  plot  were  used  for  the  equivalent  gradient 
study  plot;  however,  the  plot  identification  differed  in  each  study  unit. 
Seventy-four  plots  were  located;  36  along  transect  number  8  and  38  along 
transect  number  10. 

B.  Sub-Plot  Location 

To  facilitate  herbaceous  sampling,  two  20x50  cm  sub-plots  were 
located  within  each  plot  of  the  gradient  and  island  study  units.  The 
first  sub-plot  was  located  by  randomly  selecting  one  10°  point  in  a  180° 
interval  whose  reference  was  parallel  to  the  transect  with  0°  being 
northward.  The  longitudinal  mid-line  of  the  sub-plot  was  used  as  the 
reference  and  the  inside  edge  was  located  at  50  cm  from  the  plot  center. 
The  other  sub-plot  was  located  180°  from  the  first. 

C.  Quantitative  Species  Value 

Canopy  coverage  (%) ,  following  Daubenmire  (1959),  for  the 
herbaceous  species  was  recorded  in  each  sub-plot  and  the  average  species 
cover  of  the  two  sub-plots,  adj listed  to  the  appropriate  mid-point  of  the 
cover  class,  was  used  as  the  quantitative  species  value  for  the  plot. 


The  cover  classes  used  were: 
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Midpoint 

Range 

1.0  = 

0  - 

2% 

4.0  = 

2  - 

6% 

15  o0  = 

6  - 

24% 

37.0  = 

24  - 

50% 

63.0  = 

50  - 

76% 

85.0  - 

76  - 

94% 

97.0  = 

94  -  100% 

Species  found  within  the  plot  but  outside  the  sub-plot  were  recorded 
as  0.5%  cover  for  the  plot  while  species  found  adjacent  to  the  plot 
were  recorded  as  0.1%.  Shrub  and  tree  canopy  coverage  "was  recorded  for 
the  plot  using  the  midpoint  of  the  sub-plot  cover  classes. 

D.  Herbage  Production  Variables  Collected  in  the  Gradient  and  Island 
Study  Units 

The  green  herbage  within  both  sub-plots  was  harvested  during 
late  July  and  sorted  in  the  field  to  a  grass  and  grass-like  sample  and 
a  forb  sample.  Each  sample  was  air-dried  in  the  laboratory  and  weighed 
to  the  nearest  gram.  Annual  herbage  weight  for  each  plot  was  obtained 
by  summing  the  grass  and  grass-like  weight  and  the  forb  weight.  These 
weights  were  then  converted  to  kilograms  per  hectare. 

E.  Environmental  and  Edaphic  Variables  Collected  in  the  Gradient  Study 
unit 

Soil  moisture  and  nutrient  samples  were  collected  on  June  4, 
1971  for  each  plot.  The  elevation  of  each  plot  above  the  lake  level  was 
also  recorded.  Soil  moisture  samples  were  placed  in  moisture  tins, 
weighed,  and  oven-dried  (63°C)  for  48  hours.  The  samples  were  removed 
from  the  oven,  re-weighed,  and  the  moisture  percentage  was  calculated. 
Soil  nutrient  samples  were  sent  to  the  Alberta  Soil  and  Feed  Testing 
Laboratory  to  determine  the  pH  and  the  amount  of  nitrogen,  phosphorous, 
potassium  and  sodium  in  the  upper  15  cm  of  soil.  Soil  temperatures  were 
recorded  at  the  center  of  each  gradient  study  plot  during  the  first  week 
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of  each  of  the  following  months:  June,  July  and  August.  A  temperature 
probe  of  15  cm  in  depth  was  used  and  readings  were  made  to  the  nearest 
one-half  degree  Centigrade.  Other  variables  recorded  were  tree  and  shrub 
density  and  depth  of  organic  soil. 

III.  Data  Analysis 
A.  Ordinational  Classification 

The  ordination  technique  presented  by  Swan  et  al.  (1969)  and 
the  numerical  classification  technique  by  Carmichael  and  Sneath  (1969) 
were  chosen  for  the  present  study  because  they  both  described  the  relation¬ 
ships  between  plots  (sites)  located  in  Euclidean  space.  The  distance 
matrix  was  the  fundamental  basis  of  each  technique.  Ordination  expressed 
this  matrix  in  continuous  axes  while  numerical  classification  expressed  the 
matrix  in  discrete  clusters.  In  the  present  study  the  combination  of 
these  two  techniques  was  termed  ordinational  classification.  Ordinational 
classification  presented  the  distance  matrix  in  discrete  clusters  and 
then  ordinated  the  intercluster  distance  matrix.  Analysis  diagrams  for 
differing  approaches  to  ordinational  classification  have  been  presented 
in  Figures  4  and  5. 

1.  Gradient  and  reconnaissance  study  units. 

The  set  of  all  plots  and  species  for  each  study  unit 
comprised  a  raw  data  matrix.  This  matrix  was  analyzed  following  the 
ordinational  classification  technique  (Figure  4).  To  analyse  section  A 
of  Figure  4,  the  Euclidean  distance  program  EUCLID  (Purchase,  1970a) 
was  obtained  and  was  edited  for  a  more  efficient  data  input.  Normal 
and  transposed  raw  data  matrices  could  be  analyzed  with  EUCLID  resulting 
in  an  interplot  distance  matrix  or  an  interspecies  distance  matrix, 
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Figure  4.  Ordinational  classification  diagram  for  the  gradient  and 

reconnaissance  study  units  showing  the  related  computer  programs. 
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Figure  5.  Ordinational  classification  diagram  for  the  island 
study  unit  showing  the  related  computer  programs. 
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respectively.  The  CLUSTR  subroutine  from  Carmichael’s  TAXMAP  (1970) 
was  edited  to  use  either  distance  matrix.  This  program  classified  the 
distance  matrix  into  discrete  clusters  with  an  intercluster  distance 
matrix  being  produced  during  the  analysis  (Section  B,  Figure  4).  The 
intercluster  distance  (cluster  center  to  cluster  center)  was  determined 
from  the  sum  of  the  two  cluster  radii  plus  the  distance  between  the 
nearest  neighbors  of  the  two  clusters.  The  cluster  radius  was  defined 
as  one-half  the  maximum  interplot  distance  which  occurred  within  the 
cluster  (Carmichael,  1970).  In  the  present  study,  all  intercluster 
distances  were  used  during  cluster  ordination;  however,  only  the 
replicated  clusters  were  placed  in  the  ordination  field.  If  separ¬ 
ation  of  the  replicated  clusters  was  not  achieved,  the  position  of  the 
cluster  on  each  axis  and  the  cluster  radii  were  rescaled  to  a  percentage 
of  the  largest  axis  position  of  a  replicated  cluster. 

To  display  the  relative  proximity  of  each  cluster  in  Euclidean 
space,  three-dimensional  orthogonal  cluster  ordination  was  used.  The 
axes  of  the  orthogonal  cluster  ordination  were  defined  by  the  SDWORD 
ordination  program  (Purchase,  1970b)  which  followed  the  technique  dis¬ 
cussed  by  Swan  et  al.  (1969)  (Section  C,  Figure  4).  As  three-dimensional 
ordination  was  difficult  to  display  on  one  page,  three  two-dimensional 
ordinations  have  been  presented.  The  two-dimensional  ordination  A 
presented  the  Euclidean  space  defined  by  the  X  and  Y  axes  while  B  pre¬ 
sented  the  X  and  Z  axes  and  C  the  Y  and  Z  axes.  The  three-dimensional 
Euclidean  space  was  achieved  by  observing  the  three  two-dimensional 
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ordinations  together  as  one  three-dimensional  ordination.  The  position 
of  the  cluster  in  the  ordination  field  was  determined  by  the  cluster 
center  with  the  cluster  sphere  being  defined  by  the  cluster  radius. 

Orthogonal  plot  ordination  axes  were  also  defined  by  the  SDWORD 
program  using  the  interplot  distance  matrix  (Section  C,  Figure  4).  To 
define  each  cluster  all  plots  within  a  cluster  were  located  in  the 
ordination  field  and  a  line  was  drawn  around  the  points  delineating  the 
area  occupied  by  each  cluster.  The  display  of  the  three-dimensional 
plot  ordination  was  accomplished  by  the  same  method  as  the  three-dim¬ 
ensional  cluster  ordination.  The  plot  ordination  was  compared  to  the 
cluster  ordination  to  check  the  accuracy  of  the  cluster  analysis. 

The  clusters  of  the  reconnaissance  study  unit  have  been 
arranged  from  the  most  mesic  to  the  most  xeric  in  the  output  of  the  raw 
data  matrix  (Appendix  1)  and  in  the  one-way  analysis  of  variance  tables. 

The  clusters  of  the  gradient  study  unit  have  been  arranged  in  order  of 
increasing  elevation.  The  plots  which  have  been  presented  first  in  each 
cluster  were  at  the  center  of  the  cluster  sphere  while  the  latter  plots 
formed  the  peripheiial  or  satellite  plots  (Appendix  1).  The  clusters 
were  most  homogeneous  near  the  center,  becoming  more  heterogeneous 
as  plots  were  added  to  the  cluster.  Because  they  were  not  replicated, 
those  clusters  which  were  formed  from  single  plots  (isolated  clusters) 
have  not  been  included  in  any  analysis  or  discussion. 

The  formation  of  vegetation  associations  in  the  gradient  and 
reconnaissance  study  units  was  accomplished  by  subjectively  analyzing 
the  plot  and  cluster  ordinations  and  defining  with  the  aid  of  field  notes, 
subsets  of  clusters  which  formed  larger  groups.  The  resulting  associations 
approximately  followed  Lewis  et  al.  (1928)  and  provided  a  concise 
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description  of  the  major  vegetational  gradient. 

2.  Island  study  unit. 

The  plots  of  the  raw  data  matrix  of  the  island 
study  unit  were  normalized  so  that  the  species  with  the  highest  canopy 
coverage  in  each  plot  was  assigned  a  value  of  100%.  The  cover  values  of 
the  remaining  species  present  in  the  plot  were  adjusted  proportionately. 
This  was  done  co  remove  the  observer  error  which  had  resulted  during 
the  collection  of  canopy  coverage  data.  The  normalized  raw  data  matrix 
was  analyzed  following  the  ordinational  classification  technique  pre¬ 
sented  in  Figure  5.  The  taxometric  mapping  program  TAXMAP  (Carmichael, 
1970)  was  used  to  define  the  clusters  of  vegetation  (Section  A,  Figure  5). 
This  program  computed  the  arithmetic  distance  between  pairs  of  plots  and 
rank-ordered  all  the  distances  into  a  proximity  vector  which  was  used 
as  the  basic  matrix  during  cluster  analysis 0  This  cluster  analysis 
procedure  presented  the  "best”  set  of  clusters  from  the  raw  data  matrix 
(Carmichael,  1970).  The  intercluster  distance  matrix  was  defined  by 
TAXMAP  following  the  same  procedure  as  described  in  the  gradient  and 
reconnaissance  study  units  as  was  the  three-dimensional  orthogonal 
cluster  ordination  (Section  B,  Figure  5).  Orthogonal  plot,  ordination 
was  not  determined  for  this  study  unit  as  an  interplot  distance  matrix 
was  not  available.  The  arranging  of  the  clusters  in  the  one-way  analysis 
of  variance  tables  was  approximately  from  the  mesic  to  xeric  while  the 
isolated  clusters  were  not  included  in  any  analysis  or  discussion. 

B.  Statistical  Analysis 

Other  matrices  were  formed  from  the  plots  and  herbage  produc¬ 
tion,  environmental  and  edaphic  variables  collected  in  the  gradient 
study  unit  and  the  plots  and  herbage  production  variables  of  the 
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reconnaissance  and  island  study  units.  One-way  analysis  of  variance 
was  performed  upon  each  variable  in  the  matrices  to  determine  if 
significant  differences  (P  <  0.05)  existed  among  clusters.  Each  cluster 
was  used  as  a  vegetational  ’’treatment"  and  the  plots  within  a  cluster 
were  the  replications  within  a  treatment.  Duncan’s  New  Multiple  Range 
Test  was  applied  to  the  cluster  means  to  determine  which  were  signifi¬ 
cantly  different  (P  <  0.05).  Non-normal  (skewed)  distribution,  present 
in  some  variables,  resulted  in  a  highly  significant  homogeneity  of  vari¬ 
ance  (Chir-square)  test  when  the  analysis  of  variance  was  performed  on 
the  raw  data.  Scheff£  (1959)  reported  that  the  analysis  of  variance  on 
a  one-way  classification  could  measure  the  E-test  using  treatments  with 
unequal  variances.  He  reported  that  if  a  variance  was  largest  for  the 
treatment  with  the  greater  replication,  then  the  F-test  was  conservative 
(observed  P  -  0.04,  expected  P  =  0.05).  If,  however,  the  largest 
variance  was  associated  with  the  lesser  replicated  treatment  then  the 
F-test  was  liberal  (observed  P  =  0.09,  expected  P  =  0.05).  Steel  and 
Torrie  (1960)  discussed  this  same  point  and  concluded  that  in  biological 
data  loss  of  sensitivity  was  unimportant.  Lindquist  (1971)  reported 
that  the  F-test  was  very  elastic  and  was  able  to  tolerate  large  treat¬ 
ment  variances  without  a  loss  of  sensitivity. 

Steel  and  Torrie  (1960)  suggested  the  use  of  the  log(X+l) 
transformation  of  the  variables  to  equalize  the  variances  prior  to 
analysis.  They  also  suggested  the  use  of  Bartlett’s  test  for  the 
homogeneity  of  variances.  It  was  found  that  the  logarithmic  transform¬ 
ation  of  the  1970  forb  weight  from  the  reconnaissance  study  unit  reduced 
the  Chi-square  value  in  Bartlett’s  test  from  136.2  (raw  data)  to  29.5 
(transformed  data)  for  46  degrees  of  freedom.  The  Chi-square  value  for 
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the  transformed  variable  was  not  significant  indicating  that  the 
variances  were  homogeneous, 

Li  (1964)  reported  that  the  arithmetic  mean  from  the  logarith¬ 
mic  transformation  was  the  same  as  using  the  geometric  mean  of  the 
original  data.  He  also  reported  that  the  geometric  mean  of  a  skewed 
distribution  was  closer  to  the  median  than  the  arithmetic  mean  which 
was  affected  by  data  in  the  tail  of  the  distribution.  Steel  and  Torrie 
(1960)  reported  that  the  arithmetic  mean  may  not  be  the  most  informative 
central  value  if  the  distribution  was  skewed.  The  means  in  the  present 
study  have  been  presented  as  geometric  rather  than  arithmetic  as  the 
cluster  mean  has  been  converted  to  the  original  scale. 

Correlations  between  the  variables  in  each  study  unit  were 
determined  using  Grobben  (1970),  In  the  gradient  study  unit,  a  multiple 
regression  (Grobben,  1970)  was  determined  for  the  dependent  variable 
(annual  herbage  weight)  using  coverage,  environmental  and  edaphic  factors 
as  independent  variables,  A  multiple  regression  was  performed  on  each 
orthogonal  axis  of  the  cluster  ordination  from  the  gradient  study  unit 
to  determine  which  variables  were  Influencing  the  cluster  distribution 
along  each  axis, 

C,  Terminology  and  Species  Statistics 

A  species  found  within  a  study  unit  exhibited  a  frequency 
(Daubenmire,  1968),  a  presence  cover  (sum  of  the  cover  divided  by  the 
number  of  plots  in  which  the  species  was  present),  a  mean  cover  (sum  of 
the  cover  divided  by  the  number  of  plots  sampled)  and  a  prominence  ind ex 
(mean  cover  multiplied  by  the  square  root  of  the  frequency  (La  Hoi, 
1964)).  Dominant  species  were  defined  as  those  with  the  highest  mean 


canopy  coverage  for  the  cluster  and  the  sub-dominant  species  had  a  mean 
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canopy  coverage  which  was  at  least  60%  of  the  dominant.  The  dominant 
and  sub-dominant  species  defined  the  cluster  name  with  sub-dominance 
displayed  by  a  dash  (-)  and  co-dominance  expressed  by  a  slash  (/). 
Species  with  a  mean  canopy  coverage  of  10%  or  more  in  the  cluster  were 
defined  as  prominent s  to  distinguish  them  from  those  species  which  were 
present  with  a  low  mean  canopy  coverage  (subordinates) . 

Species  in  the  raw  data  matrix  were  arranged  by  the  computer 
in  decreasing  order  of  mean  cover  for  each  study  unit  and  an  arbitrary 
point,  approximately  1%  of  the  maximum  mean  cover,  was  used  to  delimit 
species  which  were  ’’important”  from  those  which  were  ’’unimportant". 
Inverse  ordinational  classification  of  the  gradient  study  unit  provided 
species  which  acted  as  single  species  clusters.  These  have  been  termed 
as  independent  species.  The  summation  of  the  columns  of  the  interplot, 
interspecies,  or  intercluster  distance  matrices  yielded  a  dissimilarity 
index  for  each  plot,  species  or  cluster,  respectively  (Bray  and  Curtis, 
1957).  The  plot,  species  or  cluster  with  the  largest  dissimilarity 
index  was  the  one  which  provided  the  greatest  separation  of  the  others 
in  che  ordination  field. 

IV.  Soil  Classification 
2 

Two  soil  pits  (1  m  )  were  dug  and  the  solum  was  classified 
into  the  appropriate  Order  and  notes  on  the  profiles  were  recorded.  One 
soil  pit  was  located  at  the  center  of  the  island  while  the  other  was 
located  on  the  north  slope  at  the  mesic  extremities  of  the  Salix  spp. 
Both  soil  pits  were  located  near  transect  number  10  (Figure  3).  Soil 
probings  (2.5x25  cm)  were  taken  at  each  island  study  plot  to  determine 
the  extent  of  horizon  development. 
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RESULTS 


I.  Gradient  Study  Unit 
A,  Ordinational  Classification 

The  raw  data  matrix  of  the  gradient  study  unit  was  of  the 
order  74  plots  by  64  species.  During  ordinational  classification,  56 
plots  were  arranged  into  eleven  appropriate  clusters.  The  raw  data 
matrix  has  been  presented  in  Appendix  1  showing  the  dominant  and  sub¬ 
dominant  species  which  defined  each  cluster. 

The  isolated  single  plot  clusters  had  quantitative  species 
combinations  unlike  any  other  plot.  The  species  combinations  which 
isolated  these  plots  have  been  presented  in  Appendix  2. 

Orthogonal  ordination  presented  the  relative  proximity  of  each 
cluster  in  three-dimensional  Euclidean  space.  Orthogonal  plot  ordination 
produced  large  discontinuities  among  groups  of  clusters  (Figure  6) . 

Four  associations  were  defined  from  clusters  which  were  in  close  proximity 
in  the  Euclidean  space.  These  associations  with  their  member  clusters 
were : 

(a)  The  shoreline  association; 

(i)  Cirsivm  arvense  cluster. 

(b)  The  low-moor  association; 

(i)  Calarno.gr ostis  negleota  cluster, 

(ii)  MEcotonalM  cluster, 

(iii)  Calarnagrostis  negleota  -  Corex  atherodes  cluster,  and 
(iv)  Carex  atherodes  cluster. 

(c)  The  woodland  association; 

(i)  Poa  pratensis  -  Salix  bebbiana  var.  perrostrata  cluster. 
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Figure  6.  The  three-dimensional  orthogonal  plot  ordination  o 
the  replicated  clusters  in  the  gradient  study  unit 
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(ii)  Poa  pratensis / Salix  bebbiana  var.  perrostrata  - 
Melilotus  alba  cluster, 

(iii)  Poa  pratensis  -  Populus  tremuloides  -  Populus  balsam - 
if  era  cluster,  and 

(iv)  Populus  tremuloides  -  Poa  pratensis  cluster. 

(d)  The  grassland  association; 

(i)  Poa  pratensis /Melilotus  alba  cluster, 

(ii)  Melilotus  alba  cluster,  and 
(iii)  Poa  pratensis  cluster. 

The  orthogonal  ordination  of  the  cluster  centers  in  Figure  7  produced 
cluster  proximities  in  three-dimensional  space  which  were  comparable 
to  the  orthogonal  plot  ordination  in  Figure  6. 

An  analysis  of  the  species  found  in  the  gradient  study  unit 
revealed  that  few  species  were  dominants  as  defined  by  the  following 
four  species  statistics:  frequency,  presence  cover,  mean  cover  and 
prominence  index.  The  statistics  of  the  important  species  has  been 
presented  in  Table  1»  The  species  with  the  larger  mean  cover  values 
appeared  to  be  most  influential  in  the  classification  of  the  study  unit. 
This  was  confirmed  by  inverse  ordinational  classification.  The  inverse 
ordinational  classification  isolated  ten  species  which  were  acting 
independently.  These  independent  species  were  Poa  pratensis ,  Calam- 
agros  u'ls  neglecta,  Populus  trermloides 3  Salix  bebbiana  var.  perrostrata, 
Melilotus  alkai  Populus  balsamifera3  Cirsium  arvense 3  Carex  atherodes3 
Sonehus  arvensis  and  Salix  petiolaris . 

Poa  pratensis  was  found  to  be  the  most  dissimilar  (independent) 
species  as  measured  by  the  dissimilarity  index  (Table  2) .  The  dissim¬ 
ilarity  index  was  highly  correlated  with  presence  cover,  mean  cover  and 
prominence  index  (Table  3).  Calamagrostis  negleota ,  Populus  trcmuloid.es 
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Table  1.  Statistics  defining  the  dominant 

species  in  the  gradient  study  unit. 


Species 

_ 

Poa  pratensis 

Calamagrcstis  negleota 

Populus  tremuloides 

Salix  behbiana  var.  p err o strata 

Mel-ilotus  alba 

Populus  bals  conifer  a 

Cirsiim  arvense 

3 

Carex  atherodes 

Conchas  arvensis 

Mentha  aroensis  var.  villosa 

Sa lix  pe tio laris 

Scirpus  validus 

Uvtica  qraeilis 

Ca lamagrostis  inexpansa 

Taraxacum  spp. 

Lactuca  pulchella 
Carex  praegracilis 
Astragalus  agrestis 
Juncus  balticus 
Poa  palusiris 
Scolochloa  festucacea 
Cares:  lasiocarpa 
Aster  hesperius 
Artemisia  biennis 
Hordeum  jubatum 


Freq. 

Percent 

Presence 

Covert 

Mean 

Cover 

Prominence 

Indexl 

55.4 

53.6 

29.7 

219.8 

51.4 

28.0 

14.4 

102.7 

24.3 

56.1 

13.6 

76.6 

32.4 

35.9 

11.6 

66.1 

29.7 

30,1 

8.9 

49.0 

21.6 

30.3 

6 . 6 

30.1 

59.5 

9.7 

5.7 

43.9 

44.6 

11.1 

5.0 

33.5 

62.2 

4,1 

2.5 

19.8 

35.1 

5.8 

2.0 

11.8 

13.5 

11.5 

1,6 

5.9 

35.1 

4.4 

1,5 

8,8 

28,4 

3.6 

1.0 

5,3 

36.5 

2.6 

1..0 

6.0 

52.7 

1.4 

0.7 

5.1 

21,6 

3.2 

0.7 

3.2 

10.8 

5o9 

0.6 

2.0 

8.1 

7.7 

0.6 

1.1 

32.4 

1.7 

0.6 

3.4 

36.5 

1.5 

0.6 

3 , 6 

24.3 

1.9 

0.4 

2.4 

17.6 

2.4 

0.4 

1.7 

36.5 

0.9 

0.3 

1.8 

16.2 

2.0 

0.3 

1.2 

24.3 

1.2 

0.3 

1.5 

■^liefer  to  terminology  on  page  29. 

2 

Predominantly  Poa  pratensis  with  some  Poa  cornpressa  and  Poa  glaucifolia. 
"Predominantly  Carex  atherodes  with  some  Carex  rostrata. 


■ 


'  *' 


*  > 
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Table  2.  Dissimilarity  indices  for  the  independent  species. 


Species 

Dissimilarity  Index^ 

2 

Poa.  pratensis 

53.4 

Populus  tremuloides 

37.0 

Ca lamagrostis  negtecta 

34.4 

Salix  bebbiana  var.  pevrostrata 

33,6 

Melilotus  alba 

28.7 

Populus  balscumifeva 

23,8 

Cirsiwn  arvense 

20.6 

Carex  athevodes 

18,6 

Salix  petiolaris 

12.4 

Sonchus  arvensis 

11.2 

Other  species 

<10.0 

Refer  to  terminology  on  page  30. 

2 

Predominantly  Poa  pratensis  with  some  Poa  oompressa  and 
Poa  g laucifo Ha , 

3 

'"Predominantly  Car  ex  atherodes  with  some  Cavex  vosbvata. 


«  ■ 


-43- 


Table  3.  The  relationship  between  the  dissimilarity 
index  (dependent  variable)  and  the  four 
statistics  for  the  independent  species 
in  the  gradient  study  unit. 


Species 

Statistic 

Regression 

Coefficient 

Correlation 

Coefficient 

Frequency 

0.074±0.26Q 

0.10 

Presence  Cover 

0,635±0.109** 

0.90** 

Mean  Cover 

1.516+0.130** 

0.97** 

Prominence  Index 

0.19810.028** 

0.92** 

**P  <  0.01 
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and  Salix  bebbiana  var.  perrostrata  had  equivalent  dissimilarities. 

The  ordination  of  the  independent  species  illustrated  the 
proximity  of  these  species  in  two-dimensional  space  (Figure  8) .  The 
largest  distance  was  recorded  between  Poa  pratencis  and  Calconagrostis 
negleota,  The  ecological  amplitude  of  Poa  pratensis  coincided  with  the 
higher  elevation,  lower  moisture  section  of  transect  numbers  8  and  10 
while  the  ecological  amplitude  of  Calamagrostis  negleota  corresponded 
with  the  lower  elevation,  higher  moisture  section  (Figure  9) .  The 
ecological  amplitudes  of  the  remaining  independent  species  could  be 
placed  within  the  ecological  amplitude  of  Poa  pratensis  or  Calamagrostis 
yiegleota  (Figures  10,  11).  Populus  tremuloides 3  Salto:  bebbiana  var. 
perrostrata^  Melilotus  alba s  Populus  balsamifera  and  Salix  petiolaris 
all  fell  within  the  range  of  Poa  pratensis .  Cirsiim  arvensis  and 
Carex  atherodes  recorded  the  highest  coverage  values  within  the  range 
of  Ca lamagvo s tis  yiegleota;  however,  they  were  also  present  as  a  relic 
within  the  range  of  Poa  pratensis .  Sonchus  arvensis  was  found  to  be 
present  along  most  of  transect  number  3  (Figure  10)  and  was  found  to  be 
prominent  in  the  low-moor  and  grassland  portions  of  transect  number  10 
(Figure  11) . 

B.  Statistical  Analysis 

1.  Correlation  and  regression. 

Simple  correlation  coefficients  between  variables 
calculated  from  the  data  recorded  for  each  of  the  70  plots  were  presented 
in  Table  4.  Variables  showing  significant  correlations  with  annual 
herbage  weight  were  grass  and  grass-like  weight,  forb  weight,  total 
herbage  coverage,  tree  density  and  pH.  Variables  correlated  with  grass 
and  grass-like  weight  were  annual  herbage  weight,  tree  density,  percent 
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Figure  8.  The  two-dimensional  ordination  of  the  independent  species  in  the  gradient  study  unit. 
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galamagrostis  negleota 


Figure  9.  Soil  moisture  and  elevation  relationships  for  the  two  gradient  transects. 
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Figure  10.  The  ecological  amplitude  of  each  independent 
species  along  transect  number  8. 


CODE 

SPECIES 

Popr 

Poa  pvatensis ^ 

Potr 

Populus  tvemuloides 

Sabe 

Salix  bebbiana  var.  pevvostvata 

Meal 

Melilotus  alba 

Poba 

Populus  balsamifeva 

Sape 

Salix  petiolavis 

Cane 

Ca lamagvos tis  negleota 

Ciar 

Civsium  awense 

Caat 

i  2 

Cavex  athevodes 

Soar 

Sonahus  awense 

n  -  northshore 
s  -  southshore 
t  -  highest  elevation 


Predominantly  Poa  pvatensis  with  some  Poa  compvessa  and  Poa  glauoi, folia. 
Predominantly  Cavex  athevodes  with  some  Cavex  vostvata. 
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Figure  11.  The  ecological  amplitude  of  each  independent 
species  along  transect  number  10. 


CODE 

Popr 

Potr 

Sabe 

Meal 

Poba 

Sape 

Cane 

Ciar 

Caat 

Soar 


SPECIES 

„  .  .  1 
Poa  pratens'us 

Populus  tremuloides 

Salix  bebbicma  var.  perrostrata 

Meliloius  alba 

Populus  balscanifeva 

Salix  petiolaris 

Calamagvostis  negleota 

Cirsium  arvense 

2 

Carex  atherodes 
Sonohus  arvense 


n  -  northshore 
s  -  southshora 
t  ~  highest  elevation 


Predominantly  Pea  pvatensis  with  some  Poa  eompressa  and  Poa  glaueifolia . 

2 

“'Predominantly  Carex  atherodes  with  some  Carex  vostvata . 
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moisture,  phosphorous,  sodium  and  pH,  Annual  herbage  weight  and  total 
herbage  cover  were  the  only  variables  signif icantlv  correlated  with  the 
forb  weight.  Sodium  and  percent  moisture  were  found  to  have  the  most 
significant  correlations  while  the  organic  soil  depth  was  not  signifi¬ 
cantly  correlated  with  any  variable. 

The  calculation  of  regression  of  independent  variables  (coverage, 

environmental  and  edaphic)  on  the  dependent  variable,  annual  herbage  weight, 

gave  the  equation  Y  =  -61.14  +  0.16X^  -  4.OOX2  +  13.04X2  where  Y  =  annual 

herbage  weight,  X^  =  total  herbage  cover  (forced  variable),  X^  =  tree 

density  and  X^  =  pH.  This  regression  equation  and  the  regression 

2 

coefficients  were  significant  at  the  1%  level  (R  =0.328). 

Using  only  the  12  cluster  means  as  observations,  the  correlation 
coefficients  between  the  herbage  production,  environmental  and  edaphic 
variables  were  determined  (Table  5).  "When  30  randomly  selected  pairs  of 
correlation  coefficients  were  selected  from  Tables  4  and  5,  it  was  found 
that  a  high  correlation  (r=0.711,  b=1.08,  P  <  0,01)  existed  between  the 
correlation  coefficients  of  the  two  tables.  Since  the  removal  of  the 
within-cluster  variation  did  not  significantly  change  the  correlation  be¬ 
tween  the  variables,  the  cluster  means  were  used  as  the  values  of  each 
variable. 

The  regression  analysis  using  Grobben  (1970)  for  each  orthogonal 
axis  was  derived  from  paired  observations  anxng  the  cluster  means  for  each 
variable  and  the  cluster  center  position  on  each  axis.  The  orthogonal  X- 
axis  of  the  cluster  ordination  was  defined  by  elevation,  shrub  density  and 
tree  density.  The  resulting  equation  was  Y  =  188.21  -  1,42X^  -  4.62X2  -  8.77X., 
where  Y  =  the  orthogonal  X-axis,  X^  =  elevation  (t=-5.26,  P  <  0.01),  Xn  ~ 
shrub  density  (t=-5,39,  P  <  0.01)  and  X^  =  tree  density  (t=-3,65,  P  <  0.01). 
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Table  5.  Correlation  coefficients  (r)  matrix  for  all  the  variables  in 
the  gradient  study  unit  as  derived  from  12  cluster  means. 


<« 

K 


(V 


i  e 
>  o 

0)  ‘H 

i— t  iJ 

W  cO 


<u 

>■4 

rH  3 

•H  U 
O  W 


•  rH  O. 

2f3  I 

C  to  H 


Pv 


>»  rH 
rH  *h  b 
vJ  O  oJ 

^  i*j  H 


fi)  H  ft 

S  *h  a 
^  o  3 
^  WH 


>s 

0)  *H 

G>  Id 

n  n 

H  U> 

.Q 

>> 
4-' 
3  *H 
*4  C3 

^  d 
oo  <i> 
.Q 

c 

M  H  *J 
ctj  -H  (X 
CO  O  O 

U  top 
o 

w  CO 
Q>  CO 
•H  CD 
O  C 
0)  JZ 
Or  O 
CO  tH 
P4 

rH  •  J-< 
ct3  ,0  O 

u  u  > 

U  O  O 

H  W  O 


»H  •  4J 
JO  rC 

d  u 
d  <y  i-i 
d  K  Q) 


no  ^ 

u  to 

O  *H 
pH  0) 


I 


CO  CO  O  rC 
CO  r/j  ^  GO 
Cfl  (fl  *H  *H 

u  u  »h  <y 

U  C  pS 


•H 

u 

d 

> 


o 

CO 

o 

o 

CM 

• 

• 

rH 

O 

o 

O 

co 

o 

rH 

CM 

o 

to 

«H 

o 

o 

o 


o  o 

O  CO 

O  rH 


O 

O 

o 


CO 

O 


* 

•X 

to 


o 

o 

o 


•X 

•X 

VO 

CM 

CO 


O  cm  r-* 

O  -<r  rH 

o  o  o 


* 

CO  CO 
<2*  vO 
<r  vo 

O  l’ 


rH 

Oh  to 
CO  O 


O 

o 

o 


o 

o 

o 


*x 

•X 

o 

o 

Oh 


•X 

AC 

CM 

ON 

00 


4: 

AC 

CO 

o 

00 


CM 

iO 


00 

OJ 

co 


r^- 

o 


CM 

co 

CM 


* 

o 

o 

vO 


* 

o 

to 


CM 

rH 


CO 

O 


o 

CO 

CO 


* 

CM 

CM 

vO 


o 

rH 

CO 

• 

I 


tO 

vO 


tH 

iO 


iO 


■X 

co 

CO 

vo 


O 

o 

o 


o  m 

o 

o  cm 


o 

o 

o 


vo 

C4 

<r 


00 

o 

CM 


o 

o 

o 


* 

■X 

o  r- 
O  in 
O 


AC 

IO 

m 


-X 

CO 

CO 

vo 


CO 

vo 

CO 


* 

o 

CM 

vo 


o 

00 


CM 

CO 


<r 

<r 

CO 


CT\ 

<? 

CO 


to 


vo 

oo 

CO 


o 


JO 

CO 


rO 

u 

o 

Ph 


0) 

00 

Cfl 

rO 

u 

<u 

n: 


d 

0 

c 

a 

< 


u 

a) 

> 

o 

u 

0) 

00 

d 

u 

Q> 

P3 


O 

H 


W 

w 

a) 

a 

u 

•H 

w 

a) 

*H 

a 

G) 

P- 

C/5 


PX 

QJ 


O 

CO 

a 

•H 

C 

C3 

oo 

V4 

O 


a  rx 


•H 

£0 

P 

0) 

Q 

rO 

O 

U 

,c 

Ui 


w 

a 

a) 

Q 

<i) 

a> 

M 


H 


O 

CO 

G) 

CJ 

P 

•O 


o 

CO 


o 

CO 


oo 

D 

< 


a) 

M 

P 


*r< 

O 

CO 


a 

o 


a 

> 

Q) 

rH 

w 


d 

<y 

o 

M 

H 

•H 


cn 

d 

o 

M 

O 

x: 

(X 

U) 

o 

.d 

(X 


AC 

00 

-sT 

vO 


00 

vO 

CO 


•X 

AC 

CM 

00 

00 


AC 

AC 

00 

cr> 


to 


O 

rH 

<r 

-3* 

CO 

o 

»H 

o 

Os 

O 

CM 

m 

to 

o 

00 

OJ 

00 

Oh 

o 

tH 

rH 

t — t 

o 

CM 

CM 

CO 

tH 

rH 

CM 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

O 

o 

o 

o 

i 

O 

i 

i 

o 

* 

AC 

CM 

Os 

r^. 

<3- 

CM 

<r 

00 

Cv 

o 

CO 

vO 

co 

vC 

CO 

co 

CO 

O 

CM 

<• 

MT 

CO 

CM 

O 

to 

o 

CM 

rH 

H- 

• 

a 

• 

• 

• 

• 

• 

* 

• 

» 

* 

r 

i 

I 

i 

i 

o 

t 

o 

i 

i 

o 

»o 

to 

rH 

O 

Os 

r^. 

oo 

r*- 

«o 

to 

Os 

-O 

CO 

CO 

00 

rH 

CO 

oo 

■<? 

CO 

CO 

CO 

H" 

co 

CM 

O 

rH 

CO 

rH 

o 

• 

• 

• 

e 

• 

« 

• 

o 

• 

• 

• 

o 

1 

o 

o 

o 

o 

i 

o 

i 

o 

1 

as 

o 

r^. 

H 

oo 

Os 

CO 

o 

as 

CO 

rH 

CM 

OS 

CO 

CO 

CM 

vO 

CM 

rH 

O 

CO 

tH 

tH 

<r 

CM 

O 

CM 

CM 

CM 

rH 

to 

CM 

o 

O 

o 

O 

O 

i 

O 

i 

i 

i 

O 

CM 

CM 

vO 

O 

to 

CM 

rH 

O 

CO 

CO 

MT 

<r 

CM 

to 

CM 

O 

CO 

O 

CO 

CO 

to 

to 

<r 

CO 

•sT 

CM 

tH 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

l 

1 

o 

o 

o 

o 

i 

V 

l 

i 

O 

Ac 

to 

to 

CO 

vO 

CM 

vO 

CM 

r** 

os 

<* 

rH 

HT 

uo 

CM 

CO 

tH 

rH 

o\ 

VO 

CM 

rH 

o 

to 

IO 

CM 

O 

O 

o 

CM 

CO 

CM 

1 

1 

o 

o 

O 

O 

o 

1 

i 

i 

i 

•K 

AC 

a\ 

r- 

as 

tH 

Os 

co 

H 

vO 

CM 

co 

to 

CM 

rH 

vO 

CO 

vo 

o 

00 

<T 

CM 

<r 

o 

O 

tH 

<r 

vo 

tH 

■vj* 

rH 

vO 

1 

i 

1 

o 

O 

o 

i 

i 

1 

1 

o 

vo 

cm 


vo 
c r\ 


CM 

co 

c^') 


o 

co 


O  CM 
vO  vO 
O 


O 

av 

co 


kt  cm 

VO  vo 
CO  ~sT 


O 

rH 

r*. 

Ac 

O 

AC 

co 

AC 

'O 

o 

rH 

CO 

o 

o 

ON 

r*N 

tH 

rH 

00 

<T 

o 

o 

rH 

co 

O 

O 

to 

o 

CO 

O 

to 

CM 

O 


CO 


O 

»o 

to 


00 

CM 

O 


* 

CO 

v-J- 

VO 


CM 

VO 


-3- 

CM 


W 

w 

d 

4J 

o 


o  W 

CP  P4 


o 

o 

V 

p« 

K 

■K 


- .  Y 


V 


-  - 


-  - 


« 


-54- 


The  squared  multiple  correlation  coefficient  (R  )  was  0.951  and  the 
equation  was  highly  significant  (F=51.5,  P  <  0.01).  Other  variables 
correlated  with  the  X-axis  have  been  presented  in  Table  6. 

The  orthogonal  Y-axis  was  defined  by  the  equation  Y  =  249.38  + 
48.76X^  52.59X2  -  O.78X2  where  Y  =  the  orthogonal  Y-axis,  X^  =  June 

soil  temperature  (t=7.01,  P  <  0.01),  X^  =  August  soil  temperature 
(t=-4.56,  P  <  0,01)  and  X^  =  phosphorous  (t=-.211,  P  <  0.05).  The 

2 

regression  equation  was  highly  significant  (F=19.0,  P  <  0.01)  and  R  = 
0.877.  Other  variables  correlated  with  the  Y-axis  have  been  presented  in 
Table  6, 

The  orthogonal  2-axis  was  defined  by  the  least  number  of  vari¬ 
ables,  The  resulting  equation  was  Y  =  -205.16  +  2.26X^  +  30.72X,,  where 
Y  =  the  orthogonal  2-axis,  X^  =  phosphorous  (t=4.67,  P  <  0.01)  and  X^  -  pH 
(t— 1.70),  The  R^  value  was  0.708  and  F  was  equal  to  10.9  (P  <  0.01). 

Other  variables  correlated  with  the  2-axis  have  been  presented  in  Table  6. 

2,  One-way  analysis  of  variance. 

One-way  analysis  of  variance  among  the  clusters  of  the 
gradient  study  unit  for  the  herbage  production,  environmental  and  edaphic 
variables  have  been  presented  in  Tables  7,  8  and  9,  respectively.  Signif¬ 
icant  F  values  were  found  for  each  variable  except  nitrogen  and  potassium 
confirming  that  statistical  differences  existed  among  clusters.  The 
variables  in  descending  order  of  significance  were  tree  density,  soil 
moisture  percent,  elevation,  shrub  density,  grass  and  grass-like  weight, 
sodium,  annual  herbage  weight,  forb  weight,  June  soil  temperature,  July 
soil  temperature,  species  richness,  soil  pH,  depth  of  organic  soil, 
phosphorous,  and  August  soil  temperature  (Tables  7,  8,  9). 

The  one-way  analysis  of  variance  among  associations  for  the 
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Table  6.  Correlation  coefficients  (r)  between  the  cluster  position 
on  each  orthogonal  axis  and  the  cluster  means  of  the 
herbage  production,  environmental  and  edaphic  variables 
in  the  gradient  study  unit. 


Variable 

X-axis 

Y-axis 

Z-axis 

Grass  and  grass-like  wt. 

0.57* 

-.04 

-.19 

Forb  weight 

0.06 

0.59* 

-.29 

Annual  herbage  weight 

0C45 

0.57* 

-.41 

Species  richness 

0.07 

0.24 

-.18 

June  soil  temperature 

-.14 

0.69** 

0.24 

July  soil  temperature 

0.10 

0.49 

0.07 

Aug,  soil  temperature 

0.28 

0.22 

0.29 

Soil  moisture  % 

0.81** 

-.35 

-.57* 

Elevation 

-.83** 

0.20 

0.69** 

Shrub  density 

-.59* 

0.30 

0.05 

Tree  density 

-.72** 

0.20 

0.06 

Depth  of  organic  soil 

0.62* 

-.29 

-.37 

Nitrogen 

0.30 

-.52 

0.05 

Phosphorous 

-.46 

-.26 

0.78** 

Potassium 

0.11 

-.47 

-.28 

Sodium 

0.77** 

-.39 

-.45 

pH 

0.36 

0.40 

-.60* 

*P  <  0.05 


**P  <0.01 
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Table  7.  A  one-way  analysis  of  variance  of  each  herbage  production 

variable  for  the  replicated  clusters  in  the  gradient  study  unit. 


Grass  & 

Annual 

No . 

Grass-like 

Forb 

Herbage 

Species 

of 

Weight 

Weight 

Weight 

Richness 

Cluster  Plots 

(kg/haD.M.)  (kg/haQM.) 

(kg/ha  ELM)  (no./plot) 

Ca lamagros tis  negleota 

13 

2 

2560a 

130bc 

2790ab 

9,6bcd 

"Ecotonal" 

13 

1475abc 

515ab 

2235ab 

10.0bcd 

Calamagrostis  negleota  - 
Car ex  atherodes 

2 

2400ab 

150bc 

2550ab 

6,5cd 

Car ex  atherodes 

2 

26153 

o 

o 

to 

2660ab 

6.0d 

Cirsium  aroense 

2 

595cde 

2715a 

3415ab 

16. 9a 

Poa  pratensis/ 

Melilotus  aZba 

2 

1240abcd 

2505a 

3350a 

10 ,5abc 

Poa  pratensis  -  Salix 
bebbiana  var,  perrostrata 

3 

993bcd 

2b0bc 

1255bc 

14.3ab 

Melilotus  alba 

3 

1110abca 

2590a 

3815a 

13.5ab 

Poa  pratensis /Salix 
bebbiana  var,  perrostrata 
-  Melilotus  alba. 

2 

1360abc 

2105a 

3615a 

9 . obcd 

Poa  pratensis  -  Populus 
tremuloides  - Populus 

d  e 

ho 

r  d 

prl 

balsamifera 

3 

105 

655 

7.3 

Poa  pratensis 

7 

176Qab 

620ab 

2600ab 

10.7abc 

Populus  tremuloides  - 

6.9cd 

Poa  pratensis 

4 

330e 

155 

460d 

Grand  Mean 

56 

1400±105 

3801185 

21951105 

9.711.4 

F  Value 

8.07*** 

4.89*** 

5.45*** 

3.51*** 

Haw  data  was  transformed  to  the  logarithmic  scale  and  the  means  presented 
were  geometric. 

* 

'Numbers  with  the  same  superscript  are  not  significantly  different  at 
P  <  0.05  level  (Duncan's  New  Multiple  Range  Test). 


***P  <  0.005 


-  ,’-c 

'*  .  n 


*V  ' 
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Table  8.  A  one-way  analysis  of  variance  of  each  environmental  variable 
for  the  replicated  clusters  in  the  gradient  study  unit . 1 
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herbage  production,  environmental  and  edaphic  variables  have  been 
presented  in  Tables  10,  11  and  12,  respectively.  The  woodland  assoc¬ 
iation  was  significantly  lower  in  annual  herbage  weight  than  the  low-moor, 
shoreline  and  grassland  associations  while  the  grass  and  grass-like  weight 
of  the  low-moor  and  grassland  associations  were  significantly  higher  than 
the  shoreline  and  woodland  associations  (Table  10) .  The  largest  forb 
weight  was  recorded  in  the  shoreline  association;  however,  this  association 
was  not  significantly  different  from  the  grassland  association.  The 
species  richness  of  the  shoreline  association  was  significantly  greater 
than  the  other  associations.  The  temperature  of  the  grassland  soil  was 
consistently  warmer  than  the  other  associations  although  the  woodland 
association  had  less  soil  moisture  and  occurred  at  a  high  elevation 
(Table  11) .  Soluble  soil  nutrients  were  found  to  be  the  lowest  in  the 
grassland  and  woodland  associations  while  phosphorous  was  the  highest 
(Table  12).  The.  soil  pH  was  found  to  be  significantly  lower  in  the 
woodland  association. 

II.  Reconnaissance  Study  Unit 
A.  Ordinational  Classification 

The  raw  data  matrix  of  the  reconnaissance  study  unit  was  of 
the  order  70  plots  by  73  species.  Fifty-nine  plots  were  grouped  into 
11  clusters  during  ordinational  classification  while  11  plots  were 
isolated  as  single  plot  clusters.  Arranging  the  plots  into  the  appro¬ 
priate  clusters,  the  raw  data  matrix  was  presented  in  Appendix  3  showing 
the  dominant  and  sub-dominant  species  defining  each  cluster.  The  isolated 
single  plot  clusters  contained  species  combinations  which  isolated  these 
plots  from  the  remainder.  The  dominant  and  sub-dominant  species  of  these 
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Table  10.  A  one-way  analysis  of  variance  of  each 

herbage  production  variable  for  the  four 
associations  in  the  gradient  study  unit.^ 


Association 

Grass  & 
Grass-like 
Weight 
(kg/ha  D0M.) 

Forb 

Weight 
(kg/ha  D.M.) 

Anmial 
Herbage 
Weight 
(kg/ha  D.M.) 

Species 
Richness 
(no. /plot) 

Shoreline 

b2 

598 

-  % 

2721a 

3424a 

16. 9a 

Lcw-moor 

2013a 

234° 

2510a 

9.2b 

Woodland 

Grassland 

597b 

1485a 

273bC 

1132ab 

925b 

3060a 

8.8b 

11. 3b 

F  Value 

14.10*** 

6.41*** 

10.05*** 

3.61** 

Raw  data  has  been  transformed  to  the  logarithmic  scale  and  the  means 
are  geometric 0 

Numbers  with  the  same  superscript  are  not  significantly  different  at 
P  <  0.05  level  (Duncan’s  New  Multiple  Range  Test). 


?  <  0.01 


P  <  0.005 


- 


Table  11 o  A  one-way  analysis  of  variance  of  each  environmental  variable 
for  the  four  associations  in  the  gradient  study  unit.l 
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Table  12. 


A  one-way  analysis  of  variance  of  each  edaphic  variable 
for  the  four  associations  in  the  gradient  study  unit . 1 


Association 

Organic 

Soil 

Depth  (cm) 

Nitrogen 
(kg /ha) 

\ 

Phosphorus 
(kg /ha) 

Potassium 
(kg /ha) 

Sodium 
(kg /ha) 

pH 

Shoreline 

2 

-i  “» ^ 

1  .  / 

7.8a 

b 

J  •  J 

596. 6a 

177. 9a 

8.3a 

Low -mo or 

2.0a 

4 ,5a 

5.3b 

559. 2a 

325. 3a 

7.9a 

Woodland 

0.5b 

2.8a 

21.  Sa 

517. la 

27.  lb 

7.4b 

Grassland 

1.0a 

2.2a 

29.8a 

374. 9b 

47. 3b 

8.0a 

F  Value 

4.68*** 

1.52 

6.62*** 

3.48* 

21.84** 

2.90* 

Raw  data  has  been  transformed  to  the  logarithmic  scale  and  the  means 
are  geometric. 

Numbers  with  the  same  superscript  are  not  significantly  different  at 
P  <  0.05  level  (Duncan’s  New  Multiple  Range  Test). 


*P  <  0.05 


**P  <  0.01 


***P  <  0.005 
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plots  have  been  presented  in  Appendix  4. 

Orthogonal  ordination  presented  the  relative  proximity  of  each 
cluster  in  three-dimensional  space.  Discontinuities  among  the  11 
clusters  were  revealed  by  orthogonal  ordination  of  the  plots  within 
each  cluster  (Figure  12) .  Relative  discontinuities  among  sets  of 
clusters  were  also  suggested.  Four  associations  were  defined  from 
clusters  which  occurred  in  close  proximity  (Figure  12) ,  These  associa¬ 
tions  and  their  member  clusters  were: 

(a)  The  shoreline  association; 

(i)  Scirpus  validus  -  Cirsium  arvense  -  Sonckus  arvensis 
cluster. 

(b)  The  low-moor  association; 

(i)  Scolochloa  festucacea  cluster,  and 
(ii)  Coorex  atherodes  -  Calamagrostis  neglecta  cluster, 

(c)  The  woodland  association; 

(i)  Poa  pratensis  -  Salix  bebbiana  var.  perrostrata  cluster, 

(ii)  Populus  tremuloides  -  Poa  pratensis /Salix  bebhiana 
var0  per-rostrata  cluster,  and 

(iii)  Populus  trernuloid.es  -  Poa  pratensis  cluster. 

(d)  The  grassland  association; 

(i)  Poa  pratensis  -  Astragalus  canadensis  cluster, 

(ii)  Poa  pratensis  -  Melilotus  alba  cluster, 

(iii)  Poa  pratensis  cluster,  and 
(iv)  Poa  pratensis  -  Taraxacum  spp.  cluster. 

The  orthogonal  ordination  of  the  cluster  centers  in  Figure  13 
produced  a  three-dimensional  space  which  was  comparable  to  the  orthogonal 
plot  ordination.  The  relative  proximity  of  the  clusters  was  maintained 
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but  the  size  of  the  clusters  was  larger  because  the  cluster  sphere 
was  defined  by  the  periphery  plots.  The  orthogonal  cluster  ordination 
explained  49%  of  the  variability  which  existed  within  the  intersite 
distance  matrix  while  the  orthogonal  plot  ordination  explained  41%. 

An  analysis  of  the  species  found  in  the  reconnaissance  study 
unit  revealed  that  very  few  species  were  dominants  as  defined  by  the 
four  species  statistics;  frequency,  presence  cover,  mean  cover,  and 
prominence  index  (Table  13) .  The  species  with  the  larger  mean  covers 
appeared  to  be  most  influential  in  the  classification  of  the  plots. 

B.  Statistical  Analysis 

1.  Correlation  and  regression. 

The  correlations  between  the  herbage  production 

variables  have  been  given  in  Table  14.  Since  grass  and  grass-like 

weight  plus  forb  weight  equalled  the  annual  herbage  weight,  significant; 

correlations  naturally  occurred  between  annual  herbage  weight  and  its 

components.  Other  significant  correlations  occurred  between  the  annual 

herbage  weight  for  1970  and  1971  and  the  grass  and  grass-like  weight  for 

1970  and  1971.  The  1970  standing  crop  was  significantly  correlated 

with  the  1970  and  1971  grass  and  grass-like  weight  and  the  1970  and  1971 

annual  herbage  weight.  The  1970  and  1971  forb  weights  were  not  correlated 

with  any  production  variable  except  the  1970  and  1971  annual  herbage 

weight,  respectively.  Simple  linear  regression  was  also  applied  to 

various  combinations  of  the  herbage  production  variables.  Excluding  the 

relationships  between  annual  herbage  weight  and  its  components,  only  one 

2 

significant  regression  equation  was  determined  where  r”  was  greater  than 


50%.  The  regression  equation  was  Y  =  7.02  +  0.57X  where  Y  =  1971  grass 
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Table  13.  Statistics  defining  the  dominant  species 
in  the  reconnaissance  study  unit. 


Species 

Frequency 

Percent-*- 

Presence 

Cover-*- 

Mean 

Cover- 

Prominence 

Index^ 

Poa  pratensis ^ 

58.6 

80.4 

47.1 

359.0 

*  ^ 

Calamagrostis  neglecta 

40.0 

49.6 

19.8 

124.7 

Carex  atherodes 

28.6 

50.4 

14,4 

76.3 

Populus  tremuloides 

21.4 

45.0 

9.6 

44.2 

Cirsium  arvense 

38.6 

23.1 

8.9 

55.2 

Melilotus  alba 

51.4 

17.3 

8.9 

63.2 

Taraxacum  spp. 

60.0 

14.0 

8.7 

67.0 

Salix  bebbiana  var. 
perrostrata 

24.3 

34.5 

8.4 

41.2 

Astragalus  agrestis 

40.0 

20.9 

8.3 

52.3 

Juncas  balticus 

30.0 

21.7 

6,5 

35.8 

Sonchus  arvensis 

32.9 

19.0 

6.3 

35.9 

Astragalus  canadensis 

2209 

26  o0 

5.9 

28.3 

Fragaria  virginiana 
var.  glocuca 

42.9 

13.2 

5.7 

37.6 

Epilobium  angustifoliwr 

21.4 

24.5 

5.2 

23.4 

Achillea,  millefolium 

48.6 

10.0 

4.9 

34.3 

S circus  validus 
£ 

10,0 

48.6 

4.9 

15.5 

Scolochloa  festucacea 

10.0 

43.9 

4.4 

14,1 

Solidago  spp. 

12.9 

28.0 

3.6 

13.0 

Aster  pansus 

34.3 

10.0 

3.4 

19.7 

Lactruca  pulchella 

42.9 

7.4 

3.2 

20.8 

Salix  petiolaris 

12.9 

22.8 

2.9 

10.4 

Populus  balsamifera 

20.0 

10.6 

2,1 

9.4 

Mentha  arverisis 
var.  villosa 

20.0 

9.0 

1.8 

8.1 

\efer  to  terminology  on  page  29. 

2 

Predominantly  Poa  pratensis  with  some 

P .  compressa 

and  P . 

glaucifolia. 

3 

Predominantly  C alama.gr  ostis  neglecta  with  some  C .  inearpansa. 
^Predominantly  Care: c  atherodes  with  some  C,  vostrata. 
^Predominantly  Astragalus  agrestis  with  some  A,  flexuosus . 
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and  grass-like  weight,  X  =  1970  grass  and  grass-like  weight  (t=8.6, 

2 

P  <  0.01)  and  r  =0.58.  A  significant  equation  was  also  found  for  the 

1970  to  1971  annual  herbage  weight;  however,  the  regression  equation 

2 

accounted  for  less  than  50%  of  the  total  sum  of  squares  (r  =0.20)  in 
the  analysis.  The  equation  was  Y  =  14.6  +  0.53X  where  Y  =  1971  annual 
herbage  weight  and  X  =  1970  annual  herbage  weight  (t=3.7,  P  <  0.01). 

2.  One-way  analysis  of  variance. 

Herbage  production  variables  were  collected  for  10 
of  the  clusters  during  the  first  week  of  August  in  1970  and  1971.  The 
analysis  presented  in  Table  15  revealed  that  significant  differences 
existed  for  most  variables  among  the  clusters  defined  from  canopy 
coverage  data.  Significant  differences  among  clusters  did  not  exist 
for  the  1971  annual  herbage  weight  because  the  apparent  differences 
among  clusters  were  masked  by  intra-cluster  variance  and  small  sample 
sizes.  The  forb  weight  generally  increased  from  the  mesic  to  the  xeric 
clusters 0  In  contrast,  the  grass  and  grass-like  weight  was  similar 
between  the  low-moor  and  grassland  clusters  hut  was  reduced  in  the 
woodland  clusters.  The  annual  herbage  weight  grand  means,  2480  kg/ha  D.M, 
in  1970  and  2580  kg/ha  D.M.  in  1971,  were  not  statistically  different. 

The  grass  and  grass-like  weight  grand  mean  decreased  (not  significantly) 
from  2000  kg/ha  D.M.  in  1970  to  1790  kg/ha  D.M.  in  1971  while  the  forb 
weight  grand  mean  increased  (not  significantly)  from  190  kg/ha  D.M.  in 
1970  to  250  kg/ha  D.M.  in  1971.  The  standing  crop  grand  mean  indicated 
that  litter  accounted  for  62%  of  the  6540  kg/ha  D.M.  available  forage 
in  the  study  unit  (Table  15) . 

The  annual  herbage  weights  of  the  low-moor  and  grassland 
associations  were  not  significantly  different  but  the  woodland  association 
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was  significantly  lower  than  the  former  associations  (Table  16).  It  was 
found  that  the  invasion  of  woody  species  reduced  the  mean  annual  herbage 
weight  by  1570  kg/ha  D.M.  or  61%.  The  forb  weight  was  found  to  be 
significantly  different  in  the  associations  sampled.  The  highest  forb 
weight  was  recorded  for  the  grassland  association  while  the  woodland 
was  the  next  highest  (Table  16).  The  forb  weight  of  the  low-moor  assoc¬ 
iation  was  the  lowest  with  a  mean  annual  production  of  15  kg/ha  D.M. 

The  grass  and  grass-like  weight  was  consistent  with  the  annual  herbage 
weight.  The  one  exception  to  this  generalization  occurred  in  1971 
when  the  grass  and  grass-like  weight  of  the  low-moor  association  was 
significantly  different  from  the  woodland  and  grassland  associations 
(Table  16). 


III.  Island  Study  Unit 
A.  Plot  Normalization 

Ordinational  classification  following  the  method  used  in  the 
reconnaissance  and  gradient  study  units  was  applied  to  the  raw  data 
matrix  and  13  clusters  were  defined  from  151  plots.  One  cluster  was 
large  and  very  heterogeneous  being  formed  from  plots  found  along  the 
ecotone  between  the  upland  and  lowland  zones  as  shown  in  Figure  14.  By 
analysing  the  heterogeneous  cluster  as  a  separate  unit,  12  sub-clusters 
were  formed.  Three  of  the  sub-clusters  were  found  to  be  similar  to  the 
clusters  earlier  defined  and  the  similar  clusters  and  sub-clusters  were 
analyzed  to  determine  if  the  relationship  between  weight  and  cover  was 
consistent  (Table  17).  It  was  found  that  the  total  herbage  cover  was 
significantly  lower  for  all  sub-clusters  even  though  the  annual  herbage 
weight  was  significantly  lower  in  only  the  Capex  adhepodes  sub-cluster. 
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Figure  14.  The  upland-lowland  ecotone. 
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Table  17.  Comparison  between  the  cluster  and  sub-cluster 
for  each  of  three  species  following  the  initial 
analysis  of  the  island  study  unit  raw  data  matrix. 


Item 

Annual 
Herb,  Wt. 
(kg/ha 
D.M.) 

Percent 

of 

Cluster 

Total 

Herbage 

Cover 

Percent 

of 

Cluster 

Ca tamagros tis  negteota  cluster 
Catamagvost-is  negteota  sub-cluster 

2 

3430a 

2545a 

74 

97.  la 
47. 5b 

49 

Care-x  atherodes  cluster 

Carex  atherodes  sub-cluster 

4930a 

2550b 

51 

88. 8a 
39. 8b 

45 

Oil's  ium  arvense  cluster 

Cirsium  arvense  sub-cluster 

4210a 

3600a 

86 

92. 7a 
44. 4° 

43 

Raw  data  was  transformed  to  the  logarithmic  scale  and  the  means  presented 
are  geometric. 

Numbers  with  the  same  superscript  are  not  significantly  different  at 
P  <  0.05  level  (Duncan’s  New  Multiple  Range  Test). 
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The  annual  herbage  weight  was  not  significantly  different  between  the 
clusters  and  sub-clusters  for  Calamagrostis  negleeta  or  Cirsium  arvense. 

It  was  concluded  that  the  coverage  values  had  been  underestimated 
for  plots  being  assigned  to  the  heterogeneous  cluster  and  that  plot 
normalization  may  partially  correct  this  problem.  Plot  normalization 
would  result  in  a  loss  of  coverage  information;  however,  a  better  class¬ 
ification  may  result  from  plots  which  have  been  normalized. 

B.  Ordinational  Classification 

The  ordinational  classification  of  the  data  matrix  in  the  island 
study  unit  following  plot  normalization,  resulted  in  the  formation  of 
16  clusters  into  which  154  plots  were  distributed.  The  data  matrix  was 
of  the  order  185  plots  by  70  species.  Thirty-one  plots  were  isolated 
during  the  analysis.  The  clusters  formed  as  defined  by  the  species  having 
the  largest  mean  canopy  coverage  (Table  18)  were: 

(a)  Sonchus  arvensis, 

(b)  Urtiea  gracilis, 

(c)  Cirsium  arvense  -  Seirpus  validus , 

(d)  Calamagrostis  negleeta  -  Cirsium  arvense , 

(e)  Calamagrostis  negleeta  -  Carex  atherodes , 

(f)  Calamagrostis  negleeta , 

(g)  Calamagrostis  negleeta  -  Scolochloa  festucaeea , 

(h)  Seoloekloa  festucaeea , 

(i)  Mentha  arvensis  var.  villosa , 

(j)  Salix  petiolarisy 

(k)  Salix  interior , 

(l)  Salix  hehbia.na  var,  perrostrata  // 1, 
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Table  18.  The  mean  normalized  coverage  (%)  of  the  cluster  for  the  important  species  in  the  island  study  unit 
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(m)  SaVlx  bebhiana  var.  perrostrata  //2, 

(n)  Populus  balsamifera 3 

(o)  Populus  tremuloides^  and 

(p)  Poo,  pratensis  --  Populus  tremuloides . 

Arranging  these  clusters  approximately  along  the  mesic  to  xeric  gradient, 
two  vegetational  zones  could  be  easily  defined  (Table  18) *  These  were 
the  lowland  zone  (moist)  and  the  upland  zone  (drier). 

Orthogonal  cluster  ordination  presented  in  figure  15  revealed 
that  the  clusters  formed  a  continuum  from  the  mesic  to  the  xeric  with  the 
smaller  clusters  usually  falling  among  the  larger  clusters.  Five  clusters 
representing  the  major  species  or  species  combinations  which  were  most 
replicated  have  been  presented  in  Figure  16.  These  major  (most  replicated) 
clusters  were: 

(a)  Cirsium  arvense  -  Soirpus  validus3 

(b)  C alama.gr  ostis  negleota  -  Car  ex  atherodes3 

(c)  Salix  petiolai‘is3 

(d)  Salix-  bebhiana  var„  perrostrata  #1,  and 

(e)  Poa  pvatensis  -  Populus  tremuloides . 

A  large  discontinuity  was  evident  between  the  three  largest  clusters, 

Cirsium  arvense  -  Scirpus  validus.  Calamagrostis  negleota  -  Carex  ather- 
od.es  and  Poa  prat&nsis  -  Populus  tremuloides .  These  clusters  were 

s 

analogous  to  the  shoreline,  low-moor  and  woodland  associations,  respectively 
of  the  reconnaissance  and  gradient  study  units.  The  three  largest  clusters 
illustrated  the  vegetative  polarity  which  existed  on  the  island.  It  was 
found  that  the  X  and  Z  axes  of  Figures  15  and  16  defined  the  greatest 
separation  of  these  clusters. 
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Figure  15 


-87- 


tsi 


-88- 


pi 

o 

•H 

P 

ttJ 

a 

•H 

T3 

P 

O 


P 

•H 

53 

3 


P  T) 
a)  d 

P  P 

w  co 

3 

rH  <T} 
O  PJ 
03 

r™5  i~i 
d  CO 

p:  *h 

o 

&0  0) 
O  43 
43  P 
P 

P  (3 
O  -H 

t — I  W 

fij  P 
<3  0! 
O  P 
•H  CO 

co  d 
53  r — I 
a?  o 

e 

•H  5m 

rO  O 

i  *i-5 

iy  d 

a.'  a 

5m 

43  <U 
P  4 
P 
d) 

43  <P 

H  O 


a)  <y  g) 
(3  (3  d 
d  co  d 
i — I  i—5  i — I 

a<  o.  ou 

>H  N  N 

*  >V-*  Si_; 

kN  rS  4^ 

<U  0)  0) 
43  4  4 

P  P  P 

53  C  53 

O  O  O 


ccj 

pq 

H 

to 

4 

4 

U 


CO 

y 
•XJ 
O 
Ss 
y 
43 
CO  H4 


JS 

^5 

*4 

v4 

y 


y 

§ 

§ 


i 


to 

3S 

a, 

Ss  « 

■  ~p 
y 
0) 
\-4 
C30 


•4 

Co 


y 

£ 


to  y 

*4  r-j 

+»  O 

CO  *4 


y 

y 

*4 

44 


I 

y 

CO 

S3 
y 

§ 
y 

I 

•4 

(fi  CSV 

r4  t-4  v-4 

14  y  .y  « 


co 

rH  CO 

=4=  "y 

•4 

y  y 

-54  ^4 

y  y 

Sn  s 

-{4  ® 

CO  ^ 

y  4^ 

tj  co 
y  y 

5\,T4 

y 
*  y, 
p  y 

S?  (3m 


co 

’4* 


y  -4rQ  Cs! 

?h  y  y  -54 

CJi  SHrC! 

y 

“-*.8 


y 

H 

Qm 

y 

y 


5o  O  CQ  to  Pl, 


vo 


cu 

p 

d 

t>0 

•H 

(h 


<3  4  0 


W 

O 

tJ 


r-5  CM  CO  sf  IA 


Figure  16A 


-89- 


Figure  16B 
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The  analysis  of  the  species  found  in  the  island  study  unit 
revealed  that  very  few  species  were  dominants  as  defined  by  the  four 
species  statistics  (Table  19).  Calamagvostis  negleotas  Carex  athevodes 
and  C-irsi-um  arve?lse  had  the  highest  frequency,  mean  cover  and  prominence 
index  values 0 
C.  Statistical  Analysis 

The  one-way  analysis  of  variance  for  each  herbage  production 
variable  has  been  presented  in  Table  20 .  Differences  among  clusters 
were  highly  significant  for  these  variables.  The  wooded  clusters  showed 
significantly  less  production  than  the  low-moor  clusters;  however, 
the  Poa  pratsTisis  -  Populus  tvemuloides  cluster  was  not  significantly 
different , 

The  one-way  analysis  of  variance  for  each  herbage  production 
variable  for  the  major  clusters  has  been  presented  in  Table  21,  Sig¬ 
nificant  differences  occurred  for  all  variables  except  forb  weight. 

The  differences  among  clusters  were  similar  for  the  grass  and  grass-iike 
weight  and  the  annual  herbage  weight  (Table  21).  The  upland  clusters 
were  significantly  lower  in  production  than  the  lowland  clusters.  The 
Poa  pvatensis  -  Populus  tvemuZoides  cluster  was  not  significantly 
different  from  the  lowland  clusters. 
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Table  19.  Statistics  defining  the  dominant 
species  in  the  island  study  unit. 


Species 

Frequency 

Percent-*- 

Presence 

Cover^ 

Mean 

i 

Cover 

Prominence 

Index! 

Ca lamagros tis  neglecta 

57.8 

49.7 

28.7 

218.6 

2 

Carex  atherodes 

63 . 2 

32.3 

20.4 

162.2 

Civsivm  arvense 

65.4 

27.8 

18.2 

147.0 

3 

Poa  pratensis 

36.2 

40.0 

14.5 

87.2 

Salix  bebbiana 
var.  perrostrata 

2801 

48.6 

13.7 

72.4 

Fopulus  tremu'loides 

20.5 

56.4 

11.6 

52.5 

Salix  petiolavis 

18.9 

59.6 

11.3 

49.0 

Scirpus  validus 

36.3 

21.5 

7.9 

47.9 

Sonehus  arvensis 

49.7 

13.9 

6.9 

48.9 

Mentha  arvensis 
var,  villosa 

48.6 

14,1 

6.8 

47.8 

P opu lus  ba Isamifera 

17.3 

38.5 

6 . 6 

27.7 

Urtiaa  gracilis 

3C.8 

16 . 6 

5.1 

28.4 

Juncus  balticus 

34.6 

12.9 

4.4 

26.2 

Scolochloa  festucacea 

18.4 

20.8 

3.8 

16.4 

Carex  lasiocarpa 

16,8 

16,8 

2.8 

11.5 

Poa  pains tris 

35  01 

7.3 

2,6 

15.1 

Me li lotus  alba 

15,7 

12.5 

2.0 

7.8 

Aster  hesperius 

32.4 

5.8 

1.9 

10.6 

Salix  interior 

3.8 

42.5 

1.6 

3.1 

Trifolium  repens 

2.7 

42.4 

1.1 

1.9 

Taraxacum  spp. 

29.2 

3.4 

1.0 

5.4 

Frag  aria  virginiana 
var ,  glauca 

17.8 

3.9 

0.7 

2.9 

'^■Refer  to  terminology  on  page  29. 

2 

Predominantly  Carex  atherodes  with  some  Carex  rostrata. 

3 

"Predominantly  Poa  pratensis  with  some  Poa  compressa  and  Poa  glaucifolia. 
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Table  20 e  A  one-way  analyses  of  variance  of  each  herbage  production 

variable  for  the  replicated  clusters  in  the  island  study  unit. 


Grass  &  Annual 

Grass-like  Forb  Herbage 

No.  Weight  Weight  Weight  Species 

of  (kg/ha  (kg/ha  (kg/ha  Richness 

Cluster  Plots  D.M.)  D.M.)  D.M.)  (no/plot) 


Sonchus  arvensis 

Urtica  gracilis 

Cirsium  arvense  - 
Scirpus  validus 

Calamagrostis  nelgecta  - 
Cirsium  arvense 

Calamagrostis  neglect a  - 
Car  ex  atherodes 

Calamagrostis  negleeta 

Calamagrostis  negleeta  - 
Soolochloa  festucacea 

Soo  loch  loo.  fes  tucacea 

Mentha  arvensis 
var,  villose. 

Salix  petiolaris 

Salix  interior 

Salix  l>ebbia?ia 
var.  perrostrata  #1 

Salix  hebhiana. 
var.  perrostrata  il 2 

Populus  balsamifera 

P opu lus  trernu loides 

Poo.  pratensis  - 
Populus  tremuloides 

Grand  Mean 

F  Value 


5 

785cde'2 

1865a 

3165abc 

8.1ab 

2 

1015bcde 

2520a 

4155ab 

8.5ab 

26 

1755abc 

915ab 

3710abc 

6.9ab 

2 

1225abcd 

470ab 

1845abc 

5.5b 

49 

3015ab 

145bc 

3345abc 

5.7b 

2 

2230abc 

oc 

2235abc 

6.5ab 

2 

3725a 

195bc 

3930 

7.9ab 

2 

2090abc 

2220a 

4345a 

8.5ab 

4 

2020abc 

850ab 

3005abc 

8.9ab 

16 

420def 

30bc 

525de 

8.8ab 

2 

730cde 

690ab 

1440bcd 

10. 8a 

9 

415def 

80bc 

535de 

10. 4a 

3 

305ef 

15bc 

315e 

8.3ab 

5 

d  pf 

440 

85bc 

550de 

10. 4a 

2 

175f 

203bc 

440e 

10. 5a 

23 

815ca 

365ab 

1280cde 

10. 4a 

154 

1300±185 

2651200 

18001130 

7 . 611 . 4 

13.17* * *** 

6.28*** 

17.27*** 

5.78*** 

Raw  data  was  transformed  to  the  logarithmic  scale  and  the  means  presented 
are  geometric. 

Numbers  with  the  same  superscript  are  not  significantly  different  at 
P  <  0.05  level  (Duncan’s  New  Multiple  Range  Test). 


***P  <  0.005 
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DISCUSSION 


1.  Species 

A.  The  Degree  of  Dominance  of  a  Species 

The  definition  of  the  clusters  appeared  to  be  closely  related 
to  the  independent  species  found  in  the  study  unit.  The  degree  of 
classification  was  related  to  the  number  of  independent  species  entering 
into  the  distance  equation.  The  more  species  entered,  the  finer  was  the 
resulting  classification  (sensu  Williams  and  Lambert,  1959).  To  measure 
the  individual  species  contribution  to  ordinational  classification, 
inverse  ordinational  classification  was  conducted  for  the  gradient  study 
unit  to  determine  which  species  were  acting  independently  or  individual- 
istically  (sensu  Gleason,  1926) .  It  was  found  that  those  species  with 
the  largest  presence  cover  *  mean  cover  or  prominence  index  usually  acted 
independently.  The  exception  was  Mentha  arvensts  var.  viPltosa  which 
was  grouped  with  the  non-independent  species.  The  index  of  dissimilarity 
of  the  independent  species  was  highly  correlated  with  the  presence  cover, 
mean  cover  and  prominence  index  of  the  same  species  (Table  3) .  Mean 
cover  was  the  statistic  which  correlated  the  most  closely  with  the  index 
of  dissimilarity  (r~0.93)  and  has  been  used  in  the  present  study  as  the 
measure  of  the  degree  of  dominance  of  a  species  in  the  study  unit.  The 
degree  of  dominance  was  a  measure  of  the  influence  of  an  independent 
species  in  ordinational  classification. 

The  species  having  large  mean  cover  values  appeared  to  be  most 
influential  in  the  determination  of  intersite  distances  because  they 
contributed  the  most  to  the  distance  measurement  between  sites.  This 
has  been  clearly  expressed  by  Pea  pratensis  and  Calamagrostis  negleota 
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in  the  gradient  study  unit  (Table  1) .  The  ecological  amplitude  of  Poa 
pratensis  corresponded  with  the  zone  of  vegetation  at  the  higher  elevations 
(upland  zone)  while  the  ecological  amplitude  of  Ca tamagro stis  negteota 
corresponded  to  the  zone  of  vegetation  at  the  lower  elevations  (lowland 
zone)  (Figures  9,  10,  11).  The  intersite  distance  between  a  site  in  the 
lowland  zone  and  one  in  the  upland  zone  would  be  large  because  of  the 
great  differences  in  the  quantitative  value  of  each  species,  i.e.,  the 
"presence-absence"  of  Poa  pratensis  and  Catamagrostis  negteota .  Sites 
within  either  zone  would  have  a  small  intersite  distance  because  the 
quantitative  species  values  would  be  similar.  However,  the  inclusion 
of  all  independent  species  refined  the  classification  producing  the 
"best"  set  of  clusters  inherent  in  the  data  (Carmichael,  1970). 

Ordination  of  the  independent  species  presented  the  dissimil¬ 
arity  relationships  between  these  species  (Figure  8)  which  ware  not 
evident  from  the  dissimilarity  index  (Table  2) .  The  largest  dissimilar¬ 
ity  (interspecies  distance)  was  recorded  between  Poa  pratensis  and 
Ca tamagro stis  negteota  (Figure  8).  This  dissimilarity  was  suggested 
by  the  ecological  amplitudes  of  Poa  pratensis  and  Catamagvostis  negteota 
(Figures  10s  11).  The  remaining  independent  species  were  found  to  be 
intermediate  between  Poa  pratensis  and  Catamagrostis  negteota .  The 
distribution  of  the  independent  species  in  the  ordination  field  was 
similar  to  the  cluster  distribution  along  the  X-axis  (Figures  6,  7). 

This  relationship  indicated  that  these  independent  species  were  influen¬ 
tial  in  the  orthogonal  classification  of  the  plots. 

B.  Relic  Species 

The  recently  declining  water  levels  of  Beaverhill  Lake  had 
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caused  a  rapid  succession  in  the  vegetation  on  the  island.  This 
unstabilized  state  of  the  vegetation  caused  a  problem  during  vegetational 
sampling  because  species  which  were  common  to  the  lowland  zone  were 
present  as  relics  in  the  upland  zone  (Table  18) .  Due  to  these  relic 
species,  stands  of  vegetation  were  not  easily  defined  and  were  difficult 
to  delineate.  In  the  present  study,  it  was  found  that  a  representative 
sample  of  the  vegetation  was  obtained  by  employing  a  restricted  random 
method  of  plot  location.  This  method  appeared  to  effectively  eliminate 
observer  errors  during  the  location  of  the  plot  and  resulted  in  an  equal 
chance  of  all  vegetation  in  the  study  unit  being  sampled. 

Ctrsium  avo&nse  and  Sonchus  arvensis  were  the  only  relic 
species  which  caused  distortion  in  cluster  ordination.  These  species 
recorded  the  highest  canopy  coverage  in  the  shoreline  association; 
however,  they  were  very  common  in  the  upland  zone  as  relics  (Table  18). 
This  occurrence  caused  the  clusters  of  the  shoreline  association  to  be 
placed  closer,  in  ordination,  to  the  clusters  of  the  upland  zone  than 
would  be  expected  (Figures  7,  13,  15).  This  distortion  caused  the 
clusters  of  the  low-moor  association  to  be  more  masic  on  the  X-axis 
than  the  shoreline  association  which  had  a  higher  soil  moisture  value 
(Table  11). 

II.  Island  Vegetation 

The  major  assumptions  made  in  the  present  study  were: 

(a)  The  species  composition  of  a  cluster  was  a  reliable 
criterion  for  determining  its  relative  proximity  to  other  clusters  in 
Euclidean  space. 

(b)  Species  with  similar  ecological  requirements  occurred  in 
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clusters  of  the  same  association. 

(c)  The  few  environmental  factors  affecting  the  distribution 
of  the  clusters  could  be  represented  by  axes  defining  the  orthogonal 
ordination  (sensu  Swan  et  al.,  1969). 

The  orthogonal  cluster  ordinations  presented  in  Figures  7,  13,  and  15 
for  the  gradient,  reconnaissance  and  island  study  units,  respectively, 
all  revealed  continuities  and  discontinuities  in  the  vegetational  struc¬ 
ture  of  the  island  and  satisfied  the  mentioned  assumptions.  In  general, 
the  X-axis  separated  the  lowland  clusters  from  the  upland  clusters  and 
defined  the  upland  and  lowland  zones.  It  wTould  appear  that  the  distri¬ 
bution  of  the  species  on  the  island  was  mainly  associated  with  elevation 
of  the  site  above  the  lake  level;  the  level  of  soil  moisture  was  closely 
associated,  since  the  cluster  position  on  the  X-axis  was  highly  correlated 
with  elevation  (r=-.83)  and  soil  moisture  (r=0„81)  means  for  each  cluster. 
The  elevation  and  soil  moisture  means  for  each  cluster  were  highly  cor¬ 
related  (r=--.83,  Table  5).  The  elevation  and  soil  moisture  data  from 
each  plot  ware  not  correlated  (r=-.17,  Table  4)  even  though  a  relationship 
appeared  to  exist  in  Figure  9.  Since  soil  moisture  'was  determined  only 
during  June,  its  true  relationship  with  elevation  may  not  have  been 
expressed  in  the  individual  plots.  In  the  present  study,  soil  moisture 
would  have  been  measured  three  times  had  it  not  been  for  heavy  June  rains 
which  rendered  further  determination  meaningless.  The  effect  of  the  June 
rains  was  clearly  illustrated  by  the  drop  of  0.5°C  in  the  July  grand  mean 
for  soil  temperature  compared  with  the  June  grand  mean  (Table  8). 

The  Y-axis  was  always  associated  with  a  subdivision  wTithin  the 
upland  zone.  This  was  readily  apparent  in  the  reconnaissance  study  unit 


*  < 


■V  -  - 

*. 


•  - 


w.  -  '  "  •  v  7 


■  :  s  i 


■  .  •  :jct 


*  .  • 


-*• 


1  : 

•  « 


-  . 

.  -  '  -  -r 


t 


f  ^v«  r  ■ 


> 


!  1  ■:  S 


-100- 


(Figure  13)  but  the  separation  in  the  gradient  study  unit  was  more  dramatic 

(Figure  7).  The  Salix  dominated  clusters  were  separated  from  both  the 

Populus  and  Poa  dominated  clusters  (Figure  7).  For  the  purposes  of  this 

study,  the  Salix  and  the  Populus  dominated  clusters  have  been  placed  in 

the  woodland  association  while  the  remainder  of  the  upland  clusters 

formed  the  grassland  association.  The  grassland  association  was  not 

definable  in  the  island  study  unit  because  the  Poa  dominated  plots  were 

included  in  the  Poa  pratensis  -  Populus  tvemuloides  cluster.  The  Y-axis 

in  the  island  study  unit  removed  the  Salix  dominated  clusters  from  the 

Populus  dominated  clusters.  The  regression  analysis  of  the  Y-axis  in 

the  gradient  study  unit  revealed  that  soil  temperature  at  15  cm  and  the 

2 

amount  of  available  phosphorous  were  highly  associated  (E“=Q.88)S  June 

temperature  was  the  most  significant  variable  (Table  6) . 

The  2-axis  separated  the  lowland  zone  into  the  shoreline  and 

low-moor  associations  (Figures  13,  16)  except  in  the  gradient  study  unit 

where  the  shoreline  association  was  represented  by  the  Civsiun  arvense 

cluster  (Figure  7).  The  Cirsium  covens e  cluster  appeared  to  be  a  relic 

of  an  old  shoreline;  it  was  in  transition  to  the  low-moor  association 

as  revealed  by  the  high  occurrence  of  low-moor  and  grassland  species 

(Appendix  1).  The  regression  of  environmental  and  edaphic  variables  upon 

the  Z-axis  revealed  that  the  amount  of  available  phosphorous  and  soil  pH 

2 

were  closely  associated  (R”=0.71);  however,  only  the  phosphorous  regression 
coefficient  was  significant. 

A.  The  Lowland  Zone 

The  lowland  zone  occurred  on  the  low  elevation,  high  moisture 
section  of  the  island.  The  gradient  study  unit  revealed  that  the  lowland 
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zone  had  significantly  higher  soil  moisture  and  sodium  levels  and  signif¬ 
icantly  lower  elevation  and  phosphorous  levels  than  the  upland  zone. 
Subdividing  the  lowland  zone  resulted  in  the  shoreline  and  low-moor 
associations  being  defined  by  the  Z-axis. 

1.  The  shoreline  association. 

The  shoreline  association  formed  the  outer  periphery 
of  the  island.  It  can  be  observed  as  the  dark  outer  band  in  Figure  1. 

The  shoreline  association  occupied  the  most  mesic  regions  of  the  island 
and  was  generally  dominated  by  the  moisture  tolerant  species  Cirsium 
arvense  and  Seirpus  validus,  Urtiea  gracilis  and  Sonchus  arvensis  ware 
also  very  prominent.  Seirpus  validus  was  dominant  on  the  mesic  sites 
of  the  association  and  was  analogous  to  the  reed-swamp  association  of 
Lewis  et  al.  (1928),  Cirsium  arvense,  Urviea  gracilis  and  Sonchus 
arvensis  were  found  in  the  drier  sites  of  the  association.  These 
species  were  not  found  in  the  reed-swamp  association  of  Lewis  et  al. 

(1928).  The  Cirsium  arvense  cluster  of  the  gradient  study  unit  revealed 
that  this  association  may  contain  low  quantities  of  available  phosphorous 
(Table  10)  and  an  abundance  of  moisture  (Table  9) „ 

2.  The  low-moor  association. 

The  low-moor  association  may  be  observed  in  Figure  1 
as  the  light  colored  vegetation  immediately  inside  the  shoreline  assoc- 
rLation.  The  dominant  species  forming  clusters  in  the  low-moor  association 
were  Calamagrostis  negleeta,  Car  ex  atherodes  and  Seolochloa  festueaeea. 

The  largest  area  of  this  association  was  dominated  by  Calamagrostis 
negleeta.  Lewis  et  al.  (1928)  and  Moss  (1955)  place  Carex  atherodes  and 
Seolochloa  festueaeea  in  the  reed-swamp  association;  however,  in  the 
present  study  these  species  were  closely  associated  with  Calamagrostis 
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negleeta  (Appendices  1,  3,  Table  18).  Scoloehloa  festueaoea  and  Carex 
atherodes  appeared  to  favor  the  moist  sites  of  the  low-moor  association. 
Other  prominent  species  in  this  association  were  Mentha  arvensis  var. 
V-illosa a  Cirsdum  arvense 3  Scdrpus  valddus  and  Calamagrostds  dnexpansa. 

An  example  of  the  low-moor  association  has  been  presented  in  Figure  17 
where  relic  stands  of  the  shoreline  association  may  be  observed  among 
stands  of  the  low-moor  association. 

Wroe  (1971)  found  that  the  Carex  atherodes  community  occupied 
the  depressions  which  held  -water  during  part  of  the  growing  season.  He 
noted  that  Carex  atherodes  favored  the  shallow  water  while  Sodrpus  spp. 
favored  deeper  water. 

B.  The  Upland  Zona 

The  upland  zone  occupied  the  inner  section  of  che  island  which 
had  the  higher  elevation  and  was  approximately  defined  by  the  mesic 
extremity  of  the  Saldx  spp,  (Figures  1,  14).  In  the  present  study,  this 
zone  has  been  arbitrarily  divided  into  two  association,  namely  the  wood¬ 
land  and  grassland  associations.  The  division  was  made  in  this  fashion 
to  test  the  reduction  in  herbage  production  resulting  from  the  invasion 
of  xfoody  species.  Orthogonal  cluster  ordination  of  the  reconnaissance 
study  unit  (Figure  13)  complimented  this  division  of  the  upland  zone. 
However,  the  orthogonal  cluster  ordination  of  the  gradient  study  unit 
(Figure  7)  and  the  island  study  unit  (Figure  15)  separated  the  Saldx  and 
the  Fopulus  dominated  clusters. 

1.  The  woodland  association. 

The  dominant  species  of  the  woodland  association  were 
Populus  tremuloideS;  Saldx  kehbiana  var.  perrostrata  and  Salix  petdolards 
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Figure  17.  The  low-moor  association. 
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while  Poa  pratensis,  Melilotus  alba  and  Populus  balsamifera  were  promin¬ 
ent.  An  example  of  the  Salix  dominated  clusters  has  been  given  in  Figure 
18  while  an  example  of  the  Populus  dominated  clusters  has  been  presented 
in  Figure  19.  The  canopies  of  the  woodland  clusters  ranged  from  a  rela¬ 
tively  open  canopy  as  shown  in  Figure  18  to  a  very  closed  and  dense 
canopy  as  shown  in  Figure  19.  Where  the  canopy  was  relatively  open  as 
shown  in  the  right  section  of  Figure  1,  the  herbaceous  stratum  was  well 
developed  and  was  primarily  dominated  by  Poa  pratensis  and  Melilotus  alba. 
Taraxacum  spp.,  Fragaria  virginiana  var.  glauca  and  Lactuca  pulchella 
were  prominent  forbs  in  the  Populus  dominated  clusters  while  Juncus 
balticus  was  prominent  in  the  Salix  dominated  clusters.  When  the  canopy 
was  closed  as  shown  in  the  left  section  of  Figure  1,  the  herbaceous 
stratum  was  very  suppressed.  Juncus  balticus  and  Melilotus  alba,  were 
prominent  under  the  Salix  canopy  while  Poa  pratensis  and  Taraxacum  spp. 
were  dominant  under  the  closed  Populus  canopy.  The  woodland  association 
of  the  gradient  study  unit  was  found  to  be  significantly  lower  than  the 
other  associations  in  the  depth  of  the  organic  soil,  soil  pH  and  soil 
moisture  and  higher  for  shrub  density  and  tree  density  (Tables  10,  11). 
Since  the  woodland  association  was  composed  of  shrubs  and  trees,  it  would 
be  expected  to  be  significantly  different  from  the  other  associations  in 
shrub  and  tree  densities. 

2,  The  grassland  association. 

The  grassland  association  occurred  as  openings 
between  stands  of  the  woodland  association  and  was  represented  in  only 
the  intensively  sampled  study  units;  the  gradient  and  reconnaissance 
study  units.  The  largest  grassland  stand  occurred  in  the  center  of  the 
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Figure  18. 


A  Sat'ix  dominated  stand  in  the  woodland  association. 


106- 


Figure  19.  A  Populus  dominated  stand  in  the  woodland  association. 
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island  (Figure  1) .  Poa  pratensis  was  the  dominant  grass  and  a  variety 
of  forbs  were  also  present.  Melilotus  atba  was  the  dominant  forb  while 
Astragalus  canadensis 3  Astragalus  agrestis 3  Taraxacum  spp.,  Laci^uca 
puchella  and  Epilobium  angustilfolium  were  prominent.  Clones  of 
Artemisia  lucoviciana  var.  gnaphalodes,  Thermopsis  rhombilfolia3  Cirsium 
arvense  and  Sonchus  arvensis  were  also  found  in  the  grassland  association. 
The  latter  two  species  appeared  to  be  relics  of  an  old  shoreline  assoc¬ 
iation.  A  relic  stand  of  Car'ex  atherodes  was  also  found  in  the  grassland 
association  growing  in  the  shade  of  a  large  rock. 

The  grassland  association  of  the  gradient  study  unit  wTas 
significantly  higher  in  June  and  July  soil  temperatures  than  the  other 
associations  (Table  11) .  Soil  moisture  in  the  grassland  association  was 
significantly  higher  than  in  the  woodland  association  and  significantly 
lower  than  the  shoreline  and  low-moor  associations  (Table  11) .  Of  the 
edaphic  variables  presented  in  Table  12,  the  only  significant  difference 
between  the  grassland  association  and  all  other  associations  was  in  its 
lower  quantity  of  potassium.  An  example  of  the  grassland  association  has 
been  presented  in  Figure  20;  the  area  illustrated  was  extensively  sampled 
in  the  reconnaissance  study  unit. 

III.  Herbage  Production 

It  was  readily  observable  from  Tables  10,  16  and  21  that  the 
annual  herbage  weightsof  the  shoreline,  low-moor  and  grassland  associations 
were  not  significantly  different.  In  general,  the  shoreline  association 
was  the  highest  producer  of  forbs  while  the  low-moor  association  produced 
the  highest  grass  and  grass-like  weight.  The  invasion  of  woody  species 
caused  a  marked  reduction  in  annual  herbage  weight  (Tables  10,  16,  21) 
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Figure  20. 


The  grassland  association. 
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but  little  change  in  the  environmental  and  edaphic  variables  was  observed 
(Tables  11,  12).  Therefore  the  sites  occupied  by  the  woodland  association 
have  been  assumed  to  not  differ  from  the  sites  of  the  other  associations 
in  potential  herbage  production.  Averaging  all  the  annual  herbage  weight 
means  given  in  Tables  10,  16  and  21  for  the  shoreline,  low-moor  and  grass¬ 
land  associations,  it  was  found  that  the  mean  annual  herbage  weight  was 
approximately  2820  lcg/ha  D.M.  while  the  mean  annual  herbage  weight  for 
the  woodland  association  was  approximately  905  kg/ha  D.M.  The  reduction 
of  annual  herbage  weight  associated  with  the  invasion  of  woody  species 
was  approximately  1915  kg/ha  D.M.  or  68%. 

The  standing  crop  was  also  reduced  by  the  invasion  of  woody 
species  (Table  16).  The  reduction  observed  was  2765  kg/ha  D.M.  or  37%. 
Litter  accounted  for  4000  kg/ha  D.M.  or  62%  of  the  potentially  available 
forage  of  6540  kg/ha  D.M.  found  in  the  reconnaissance  study  unit  (Table 
15) . 

"Wroe  (1971)  found  that  communities  in  Central  Alberta  comparable 
to  the  low-moor,  woodland  and  grassland  associations  had  an  annual 
herbage  production  of  5323,  520  and  1600  kg/ha  D.M.,  respectively. 

Hilton  (1970)  using  similar  Central  Alberta  community-types  found  approx¬ 
imately  3800,  500  and  1100  kg/ha  D.M.,  respectively.  The  woodland  and 
grassland  associations  in  the  present  study  appeared  to  be  comparable 
with  the  results  obtained  by  Wroe  (1971)  and  Hilton  (1970)  (Tables  10,  16, 
21).  The  low-moor  association  of  the  present  study  however,  appeared  to 
be  less  productive  than  that  found  by  Wroe  (1971)  and  Hilton  (1970).  No 
communities  comparable  to  the  shoreline  association  were  studied  by 
Uroe  (1971)  or  Hilton  (1970). 
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The  annual  herbage  weight  for  each  plot  was  found  to  be  signif¬ 
icantly  correlated  with  soil  pH  (r=0.47),  tree  density  (r=-.38),  and 
total  herbage  cover  (r=0.26)  of  that  plot  (Table  4)  while  the  annual 
herbage  weight  cluster  mean  was  correlated  with  soil  pH  (r=0.64)  and 
total  herbage  cover  (r=0.56)  means  (Table  5).  Since  soil  pH  was  reduced 
following  the  invasion  of  Populus  tvemuloides  (Dormaar  and  Lutwick,  1966), 
the  reduction  of  the  annual  herbage  weight  appeared  to  be  directly 
attributable  to  the  density  of  the  trees. 

The  plot  grass  and  grass-like  weight  was  significantly  corr¬ 
elated  with  tree  density  (r=-.44),  sodium  level  (r-0,42),  soil  pH  (r=0.41) 
and  soil  moisture  (r=Q.4Q)  (Table  4).  The  cluster  grass  and  grass-lika 
weight  mean  was  only  significantly  correlated  with  the  sodium  level 
(r=0.65)  and  elevation  (r=-.63).  Since  the  grass  and  grass-like  species 
weight  was  a  component  of  the  annual  herbage  weight,  it  also  appeared  to 
be  influenced  by  the  tree  density. 

The  environmental  and  edaphic  variables  which  were  significantly 
correlated  with  the  tree  density  were  soil  moisture  (r=-.49),  soil  pH 
(r=-o40),  August  soil  temperature  (r=-.36),  sodium  (r=-.31)  and  July  soil 
temperature  (r=-.27)  (Table  4).  The  tree  density  was  also  significantly 
correlated  with  the  X-axis  (r=-.72)  (Table  6).  Therefore  it  would  appear 
from  this  data  that  the  invasion  of  the  woody  species  as  expressed 

by  tree  density  had  a  significant  influence  upon  the  classification  of 
vegetation  and  the  reduction  of  potential  herbage  production  in  the 
woodland  association. 
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IV.  Soils  of  the  Study  Area 

Soils  on  the  island  were  found  to  be  of  the  Regosolic  Order. 
Horizon  development  was  not  detected  in  the  two  soil  pits  dug  and  prob¬ 
ings  at  each  plot  in  the  island  study  unit  failed  to  show  any  substantial 
development.  The  texture  of  the  parent  material  consisted  of  three 
types;  fine  sand,  coarse  sand  and  fine  gravel,  and  clay  (Figure  21).  The 
upland  soil  was  classified  as  Orthic  Regosol  having  a  0-5  cm  LH  horizon 
above  a  C  horizon.  The  color  of  the  C  horizon  in  the  till  was  dark 
greyish  brown  (10YR4/2m).  The  lowland  soil  was  also  classified  as  Orthic 
Regosol;  however,  Cumulic  Regosol  and  Gleyed  Regosol  profiles  were  found. 
It  was  assumed  that  the  vegetation  has  not  had  a  profound  effect  upon 
soil  development  and  the  vegetation  and  soil  have  not  reached  an  equil¬ 
ibrium  since  the  recent  exposure  of  the  island  above  the  lake  level  some 
40  years  ago. 

Species  which  were  more  common  in  the  eastern  section  of  the 
island  were  Populus  tremuloides ,  Poa  pratensis3  Trifolium  repens3  Cirsium 
arvense  and  Sonohus  arvensis.  Populus  balsamifera 3  Salix  bebbiana  var. 
perro strata ,  Salix  petiolarisj  Calamagrostis  negleota  and  Carex  atherodes 
were  most  common  in  the  western  section.  Melilotus  alba  was  most  common 
in  the  center  of  the  island.  Stands  of  the  woodland  association  were 
much  more  dense  In  the  western  section  than  the  east  (Figure  1).  Herbage 
production  was  much  reduced  in  all  woodland  stands  of  the  western  section 
while  Poa  pratensis  could  be  found  forming  a  prominent  stratum  under  the 
relatively  open  woodland  canopy  in  the  eastern  section.  This  difference 
in  woodland  stands  appeared  to  be  related  to  the  texture  of  the  parent 
material.  The  coarse  sand  and  fine  gravel  would  likely  provide  a 


' 


*  A 


. 


,v 


Figure  21.  The  texture  of  the  parent  material  on  the  study  island. 
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continuous  water  table  at  about  60  cm.  Nutrient  status  would  probably 
be  the  limiting  factor  for  growth  in  the  eastern  section.  The  western 
section  was  underlain  by  high  sodium  clay  and  water  movement  from  the 
lake  was  probably  non-existent.  Therefore  the  limiting  factor  for  growth 
in  the  western  section  would  probably  be  moisture. 

The  division  in  the  woodland  vegetation  between  the  eastern 
and  western  sections  did  not  coincide  with  the  division  between  the 
coarse  sand  and  fine  gravel  and  the  clay  sections.  The  vegetation 
division  occurred  some  200  m  to  the  west  of  the  soil  division.  This 
200  m  section  was  composed  of  the  eastern  woodland  type  surrounding  a 
central  grassland  in  which  Melilotus  alba  and  Poa  pratensis  were 
dominants  (Figure  1) . 

Soil  moisture  analysis  of  the  upper  15  cm  revealed  the  PopulliS 
tremuloides  -  Poa  pratensis  and  the  Poa  pratensis  -  Populus  tv emu loides  - 
Populus  balsamifera  clusters  were  on  drier  soil  than  the  other  clusters 
(Table  8).  Transect  number  8  was  drier  than  transect  number  10  in  the 
middle  of  the  island  (Figures  10,  11).  Transect  number  8  passed  through 
stands  of  the  Populus  tremuloides  -  Poa  pratensis  and  Poa  pratensis  - 
Populus  trenruloides  -  Populus  balsamifera  clusters  while  transect  number 
10  passed  through  stands  of  the  Poa  pratensis  and  Poa  pratensis  -  Melilotus 
alba  clusters.  Soil  moistures  from  the  few  plots  of  transect  number  8 
which  were  found  in  the  coarse  sand  and  fine  gravel  textured  parent 
material  revealed  that  the  upper  15  cm  were  drier  than  the  clay  textured 
parent  material  immediately  adjacent.  A  stand  of  the  Populus  tremuloides  - 
Poa  pratensis  covered  both  types  of  soil  textures.  This  drier  condition 
was  probably  accounted  for  by  the  excessive  drainage  of  the  coarse 
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textured  parent  material.  Although  soil  texture  appeared  to  be  subtly 
influencing  the  distribution  of  the  vegetation,  further  work  would  need 
to  be  done  to  determine  cause  and  effect  relationships  which  exist 
between  soil  texture  and  vegetation. 

V.  Ordinational  Classification 

A.  x\dvantages 

Daubenmire  (1952)  suggested  that  the  study  of  the  vegetative 
structure  was  a  practical  means  of  recognizing  ecosystems.  He  also 
suggested  that  it  was  more  efficient  to  study  vegetation  using  conscious 
stand  selection.  Greg-Smith  (1964)  also  reported  that  stand  selection 
was  profitable  but  cautioned  that  conscious  stand  selection  may  present 
discontinuities  in  the  data  which  in  reality  did  not  exist.  Bray  and 
Curtis  (1957)  suggested  that  vegetation-environment  relationships  should 
be  derived  from  a  direct  analysis  of  the  quantitative  vegetational  data 
rather  than  from  stands  which  had  been  previously  classified.  In  their 
quest  to  find  the  relationships  between  stands,  Bray  and  Curtis  (1957) 
ignored  the  basic  structure  of  the  vegetation,  i.e.,  the  community. 

This  information  had  to  be  extrapolated  from  the  ordination. 

In  ordinational  classification,  an  attempt  has  been  made  to 
present  the  botanical  information  of  community  structure  in  the  relative 
continuity  or  discontinuity  of  the  vegetation  inherent  in  the  study  unit. 
Ordinational  classification  appeared  to  combine  the  community  structure 
of  Daubenmire  (1952)  and  the  stand  relationships  of  Bray  and  Curtis  (1957). 
This  approach  delimited  the  clusters  (communities)  and  displayed  the 
relationships  among  the  clusters. 

In  the  present  study,  the  conscious  stand  selection  was  attempted 
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in  the  reconnaissance  study  unit;  however,  the  successional  stage  of  the 
vegetation  on  the  island  rendered  stand  selection  impractical.  A  restricted 
random  method  of  plot  location  resulted  in  all  vegetation  having  an  equal 
opportunity  of  being  sampled.  This  has  been  clearly  illustrated  in  the 
clusters  of  the  island  study  unit.  The  major  clusters  were  represent¬ 
ative  of  the  predominant  vegetational  groupings  and  appeared  to  be 
analogous  to  those  communities  which  would  have  been  formed  through 
conscious  stand  selection,  i.e.,  following  the  method  of  Daubenmire 
(1952),  On  the  other  hand,  the  minor  clusters  would  probably  have  been 
ignored  during  conscious  stand  selection  as  they  appear  to  be  represent¬ 
ative  of  the  ecotones  between  the  major  clusters.  The  Sonahus  arvensis 3 
Mentha  ar-vensis  var,  vittosa  and  Populus  halsamifeva  clusters  seemed  to 
represent  the  adaptation  of  a  species  to  a  localized  habitat  (Table  13). 

The  stands  of  these  three  clusters  were  found  widely  dispersed  on  the 
island  and  represented  a  small  percentage  of  the  vegetation  (Table  20). 
The  other  minor  clusters  presented  an  example  of  the  ecotones  which 
existed  between  the  major  clusters  (Table  18).  By  including  the  minor 
clusters  in  the  orthogonal  cluster  ordination,  the  vegetational  continuum 
inherent  in  the  data  was  expressed  (Figure  15).  Rejection  of  these  minor 
clusters  produced  large  discontinuities  between  the  major  clusters 
(Figure  16),  Grdinational  classification  presented  a  more  complete 
expression  of  the  vegetational  structure  of  the  island  study  unit  than 
either  plot  ordination  or  conscious  stand  selection  because  ordinational 
classification: 

(a)  defined  subsets  of  similar  plots  (clusters)  (Appendices  1,  3), 

(b)  presented  the  community  structure  inherent  in  the  data 
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(Table  18),  and 

(c)  maintained  the  relative  positions  of  the  cluster  during 
cluster  ordination  (Figures  6,  7,  12,  13). 

B.  Problems 

A  problem  similar  to  that  encountered  by  Swan  et  al.  (1969) 
was  experienced  in  the  present  study.  If  no  species  clearly  dominated 
the  plot,  this  plot  could  not  be  properly  classified.  Plots  which  had 
no  dominant  species  tended  to  form  heterogeneous  clusters,  i.e.,  the 
"Eco tonal”  cluster  in  Appendix  1.  This  heterogeneous,  low-cover  cluster 
resulted  from  relatively  short  intersite  distances  between  plots  in  the 
cluster.  The  short  distances  were  not  a  result  of  plot  similarities 
but  were  a  result  of  relatively  low-cover  values  of  species  which  dom¬ 
inated  the  plots.  The  majority  of  these  "low-cover”  plots  were  found 
along  the  ecotone  between  the  upland  and  lowland  zones.  An  example  of 
this  ecotone  has  been  given  in  Figure  14.  The  canopy  coverage  of  the 
dominant  herbaceous  species  usually  occupied  less  than  50%  of  the  plot 
area  while  the  remainder  was  litter  or  bare  ground.  The  reduction  in 
the  canopy  coverage  was  probably  caused  by  the  competition  for  soil 
nutrients  and  moisture  by  the  roots  of  the  adjacent  Sal'lx  plants.  It 
was  observed  that  the  root  zone,  as  measured  by  the  reduction  In  canopy 
coverage  of  herbaceous  species,  was  much  larger  than  the  coverage  of 
the  Sal'lx  canopy  (Figure  19)  .  Plots  falling  within  the  Sal'lx  root  zone 
but  outside  the  canopy  were  highly  variable;  species  composition  changed 
erratically  and  annual  herbage  weight  ranged  from  250  lcg/ha  D.M.  to 
3200  kg/ha  D.M. 

In  the  present  study,  these  "low-cover”  plots  distorted 
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classif ication.  Sines  there  were  only  two  transects  in  the  gradient 
study  unit,  few  T,low-covern  plots  were  encountered  and  these  have  been 
left  as  a  heterogeneous  cluster  ("Ecotonal"  cluster,  Appendix  1).  In 
the  island  study  unit,  plot  normalization  eliminated  the  problem  of 
’'low-cover"  plots  because  it  scaled  the  cover  values  of  the  dominant 
species  in  each  plot  to  100  and  proportionally  adjusted  the  remaining 
cover  values.  This  method  equalized  the  cover  values  of  the  dominant 
species  in  each  plot  so  that  each  plot  entered  ordinational  classific¬ 
ation  with  a  common  base.  Plot  normalization  did  not  appear  to  change 
the  basic  structure  of  the  vegetation  (Table  18);  however,  the  basic 
canopy  coverage  data  was  lost. 

VI.  Relevance  to  Range  Management 
Ordinational  classification  has  provided  a  method  for  delin¬ 
eating  groupings  of  vegetation  which  can  be  subjected  to  statistical 
analysis.  Since  herbage  production  was  the  primary  interest  in  this 
study,  significant  differences  have  been  shown  between  clusters  for  the 
herbage  production  variables  as  well  as  between  associations.  Range 
management  decisions  can  be  based  upon  these  differences. 
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SUMMARY  AND  CONCLUSIONS 


The  synthesis  of  ordination  and  cluster  analysis  in  the  present 
study  has  been  termed  ordinational  classification.  Ordinational  class¬ 
ification  presented  the  community  structure  of  the  vegetation  and  the 
relationships  among  groupings  of  vegetation.  Two  methods  for  determining 
the  intersite  distance  have  been  presented.  These  were  the  Euclidean 
and  Arithmetic  distances.  Problems  encountered  from  low-cover  dominant 
species  have  been  overcome  in  the  present  study  by  using  plot  normaliz¬ 
ation.  It  was  more  advantageous  to  sample  the  successional  vegetation 
using  a  method  of  restricted  random  plot  selection  rather  than  conscious 
stand  selection. 

It  was  found  that  presence  cover,  mean  cover  or  prominence 
index  were  reliable  indicators  of  the  influence  of  an  independent  species 
in  ordinational  classification.  Mean  cover  was  most  closely  correlated 
with  the  dissimilarity  index  of  the  independent  species.  Poa  pratensis 
and  Calccnagvost-is  negleota  were  the  two  most  dissimilar  species.  Poa 
pratensis  defined  the  upland  zone  and  Calamagvost-Ls  negleota  defined  the 
lowland  zone.  The  other  independent  species  were  found  to  fall  Into 
either  of  the  two  zones. 

The  orthogonal  axes  defined  during  ordination  'ware  found  to  be 
closely  associated  with  environmental  and  edaphic  factors.  The  first 
(X)  orthogonal  axis  was  significantly  correlated  with  elevation,  soil 
moisture,  sodium,  tree  density,  depth  of  organic  soil,  shrub  density 
and  grass  and  grass-like  weight.  During  regression  analysis,  95%  of 
the  variation  in  the  X-axis  was  accounted  for  by  elevation,  shrub  density 
and  tree  density.  The  X-axis  defined  the  upland  and  lowland  zones. 
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The  second  (Y)  orthogonal  axis  was  significantly  correlated 
with  June  soil  temperature,  forb  weight  and  annual  herbage  weight. 
Regression  analysis  of  the  Y-axis  included  June  soil  temperature,  August 
soil  temperature  and  phosphorous  as  the  significant  factors.  These 
factors  accounted  for  80%  of  the  variation  in  the  Y-axis. 

The  third  (Z)  orthogonal  axis  was  significantly  correlated 
with  phosphorous,  elevation,  soil  pH  and  soil  moisture.  The  regression 
analysis  which  accounted  for  71%  of  the  variation  of  the  Z-axis  included 
phosphorous  and  soil  pH  as  defining  factors.  Phosphorous  was  the  only 
significant  factor. 

Three  vegetation  associations  were  defined.  These  were  the 
shoreline,  low-moor  and  woodland  associations.  A  fourth  association, 
the  grassland  association,  was  also  defined  when  the  study  unit  was 
intensively  sampled.  The  shoreline  and  low-moor  associations  were  found 
within  the  lowland  zone.  These  associations  were  defined  by  the  Z-axis. 
The  woodland  and  grassland  associations  were  found  in  the  upland  zone. 
This  division  was  suggested  by  the  Y-axis;  however,  an  arbitrary  separ¬ 
ation  was  made  in  the  present  study  to  test  the  reduction  of  the  poten¬ 
tial  herbage  production  following  the  invasion  of  woody  species. 

The  annual  herbage  productiors  of  the  shoreline,  low-moor 
and  grassland  associations  were  not  significantly  different.  The  wood¬ 
land  association  was  significantly  less  productive.  The  invasion  of 
woody  species  reduced  the  annual  herbage  production  by  approximately  68%. 

The  soils  in  the  study  area  were  of  the  Regosolic  Order.  The 
vegetation  has  not  had  a  profound  effect  upon  soil  development  and  an 
equilibrium  between  vegetation  and  soil  has  not  been  reached. 
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Appendix  1.  Canopy  coverage  (%)  of  important  species  defining  the  replicated  clusters  in  the  gradient  study  unit. 
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Appendix  2.  Dominant  and  subdeminant  species  defining  the  single 
plot  clusters  in  the  gradient  study  unit. 
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Appendix  3.  Canopy  coverage  (%)  of  important  species  defining  the  replicated 
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Appendix  3.  Canopy  coverage  (%)  of  important  species  defining  the  replicated  clusters  in  the  reconnaissance  study  unit. 
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Appendix  4.  Dominant  and  subdominant  species  defining  the  single 
plor.  clusters  in  the  reconnaissance  study  unit. 
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Appendix  5.  Botanical  and  common  names  of  plant  species 


Achillea  millefolium 

Common  yarrow 

Artemisia  biennis 

Biennial  wormwood 

Artemisia  ludoviciana  var.  gnaphalodes 

Prairie  sage 

Aster  hesperius 

Western  willow  aster 

Aster  pansus 

Many-flowered  aster 

Astragalus  agrestis 

Purple  milk-vetch 

Astragalus  canadensis 

Canadian  milk-vetch 

Astragalus  flexuosus 

Slender  milk-vetch 

Bechnannia  syzigaehne 

Slough  grass 

Calamagrostis  amerioana 

B.ead  grass 

Calamagrostis  canadensis 

Marsh  reed  grass 

Ca lamagros tis  inexpansa 

Northern  reed  grass 

Calamagrostis  n&glecta 

Narrow  reed  grass 

Carex  atherodes 

Awned  sedge 

Car ex  lasiocarpa 

Woolly  sedge 

Carex  praegraci lis 

Graceful  sedge 

Carex  rostrata 

Beaked  sedge 

Cirsium  arvense 

Canada  thistle 

Epilobium  angus Ufa  Hum 

Fireweed 

Festuca  scabrella 

Rough  fescue 

Frag  aria  virginiana  var.  glauca 

Wild  strawberry 

Glyceria  grandis 

Manna  grass 

Hordeum  jubatum 

Wild  barley 

Juricus  balticus 

Baltic  rush 

Lactuca  pulchella 

Blue  lettuce 
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Appendix  5. 

Continued. 

Melilotus  alba 

White  sweet  clover 

Mentha  arvensis  var,  villosa 

Wild  mint 

Phalaris  arundinacea 

Reed  canary  grass 

Phragmites  communis 

Giant  reed  grass 

Poa  compressa 

Canada  bluegrass 

Poa  glaucifolia 

Glaucous  bluegrass 

Poa  palustris 

Fowl  bluegrass 

Poa  pratensis 

Kentucky  bluegrass 

F o lyaonum  avicu lave 

Common  kno tweed 

Populus  bats  coni  f era 

Balsam  poplar 

Populus  tremuloides 

Aspen  poplar 

Rumex  spp. 

Sorrel 

Salix  bebbiana  var,  perrostrata 

Beaked  willow 

Salix  interior 

Sandbar  willow 

Salix  petiolaris 

Pussy  willow 

Scirpus  validus 

Great  common  bullrush 

Scclockloa  festucacea 

Whitetop 

Sonchus  arvensis 

Sow  thistle 

Taraxacum,  officinale 

Common  dandelion 

Thermopsis  rhombi folia 

Golden  bean 

Trifolium  repens 

Alsike  clover 

Typha  latifolia 

Common  cattail 

Urtica  gracilis 

Common  nettle 
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